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1.0  TOTAL  INTERNAL  PARTITION  FUNCTIONS 


The  knowledge  of  the  total  internal  partition  sum  (TIPS)  is 
useful  for  a  number  of  reasons.  For  a  system  in  thermodynamic 
equilibrium,  the  states  obey  a  Boltzmann  distribution.  Hence,  the 
number  density  of  molecules  in  a  particular  state  at  a  given 
temperature  can  be  obtained  from  the  total  internal  partition 
function.  Relationships  between  the  intensity  of  a  spectral  line,  the 
transition  moment  squared,  Einstein  coefficients,  or  oscillator 
strengths  require  the  partition  sum  be  known.1  >2  if  the  intensity  of  a 
spectral  line  is  known,  the  intensity  at  other  temperatures  can  be 
obtained  knowing  the  partition  sum  at  both  temperatures.  This  has 
important  consequences  when  studying  systems  that  are  not 
isothermal,  e.g.  the  terrestrial  atmosphere.  Having  a  database  of 
spectral  lines  of  all  important  atmospheric  species3’4  and  knowing 
the  partition  sum  at  atmospheric  temperatures,  atmospheric  spectra 
can  be  inverted  to  obtain  concentration  profiles.  Note,  in  the 
literature,  the  term  partition  function  is  synonymous  with  the  term 
partition  sum. 

The  calculations  of  the  total  internal  partition  sums  for  the  list 
of  species  taken  from  the  HITRAN  database3  were  made  under  this 
project.  However,  the  determination  of  the  TIPS  is  not  enough.  In 
order  to  be  useful  for  reduction  of  spectra  a  rapid  yet  accurate 
means  of  determining  the  partition  sum  is  needed. 

The  total  internal  partition  sum  is  defined  as  a  direct  sum  over 
all  states  of  the  molecule  of  the  factor  exp(-Ej/kT).  This  can  be 
simplified  by  accounting  for  degenerate  states  by  the  formula 


1 


E:/kT 


0) 


Q(T)  -  Sd,d,e- 

all  states 


where  dj  accounts  for  degeneracies  which  are  state  dependent.  An 
example  of  such  degeneracies  can  be  found  using  H2160.  Because  of 
the  C2V  symmetry  of  the  molecule  there  is  a  coupling  between 
nuclear  states  of  the  hydrogens  and  rotational  states.  In  the  ground 
vibrational  state,  the  states  of  odd  and  even  Ka+Kc  have  degeneracy 
factors  of  3  and  1,  respectively.  Thus,  the  degeneracy  factor  depends 
on  the  rotational  state.  The  factor  dj  accounts  for  state  independent 
degeneracy  factors.  An  example  of  such  factors  are  the  nuclear 
degeneracy  factor  that  do  not  couple  with  rotational  motion  (e.g.  non- 
symmetrical  molecule  such  as  CO).  This  gives  rise  to  a  degeneracy 
factor  given  by  the  product  of  (2Ij+l),  where  Ij  is  the  nuclear  spin  of 
nuclei  j.  Thus,  for  12C160  the  state  independent  degeneracy  factor  is 
(2  l/2+l)(2  0+1  )=2.  The  other  terms  in  the  equation  are:  E;  is  the 
total  energy  (electronic,  vibrational,  rotational,  and  any  other 
quantized  motion)  of  state  i,  k  is  the  Boltzmann  constant,  and  T  is  the 
temperature  in  K. 

In  order  to  calculate  the  partition  sum  the  degeneracy  factors 
and  energy  levels  of  the  species  in  question  are  needed.  As  the 
energy  of  the  state  increases,  the  exponential  factor  approaches  zero 
and  the  sum  can  be  truncated  with  no  loss  of  accuracy.  Where  this 
occurs  depends  on  the  temperature  since  kT  divides  the  energy  of 
the  state.  Accurate  calculations  of  the  partition  sum  can  only  be 
made  by  summing  over  all  energy  levels  until  the  exponential  factor 
no  longer  contributes  to  the  sum.  This  can  be  determined  by 


2 


considering  Q(T)  versus  energy  and  noting  where  Q(T)  does  not 
change  as  additional  energies  are  summed  over.  At  this  point  the 
calculation  of  the  partition  sum  is  said  to  have  converged. 

There  are  many  conditions  for  which  energy  levels  are  not 
known  to  high  enough  levels  to  calculate  accurate  partition  sums. 
For  these  cases  approximations  must  be  made.  One,  called  the 
product  approximation,  can  be  made  when  many  energy  levels  are 
available  for  the  ground  rotational  state  but  those  for  excited 
vibrational  states  are  not  known.  The  assumption  amounts  to 
expressing  the  energy  for  the  excited  vibrational  states  as  the  sum  of 
the  vibrational  and  rotational  energies,  Evj  =  Ev+Ej.  Thus,  Eq.  (1)  can 

be  written 


Q(T)  = 


X  dvib  e'Evii,/kT  £drole-E"-'kT 

vibrational  rotational 

slates  states 


Q  vib  Q rot 


(2) 


where  dvib  and  drm  are  the  degeneracies,  and  Evib  and  Erot  are  the 
energies  of  the  vibrational  and  rotational  states,  respectively.  When 
the  energies  of  the  rotational  states  are  not  known  to  high  enough 
levels  for  the  temperatures  in  question,  there  are  a  number  of 
analytical  formulas  of  varying  degrees  of  accuracy  that  can  be  used. 
These  range  from  crude  classical  formulas5  to  sophisticated  analytical 
developments.6" 1 1  For  the  species  under  question  the  following 
approach  was  taken.  When  all  energy  levels  are  available  such  that 
the  partition  sums  are  converged,  the  direct  sum  is  used.  When  the 
energy  levels  of  the  upper  vibrational  states  are  not  known  but  the 
rotational  levels  for  the  ground  vibrational  state  are,  the  product 


3 


approximation  is  used.  Finally,  when  there  are  a  significant  number 
of  energy  levels  missing  for  the  temperature  under  consideration, 
the  best  analytical  formulas  are  used. 

Once  the  partition  functions  are  calculated,  it  is  useful  to 
organize  them  in  a  form  that  is  convenient  for  rapid  accurate  recall. 
While  interpolation  methods  can  accomplish  this,  they  require  the 
storage  of  many  data  and  can  be  quite  slow.  A  polynomial 
expression  in  T  has  been  found  to  reproduce  the  partition  sums  to 
high  accuracy  for  most  of  the  species  on  the  HITRAN  database  and 
provides  a  rapid  means  of  retrieving  the  partition  sum.  The 
expression  is  3r(1  order  in  temperature  and  requires  the  storage  of  4 
coefficients  for  a  given  temperature  range,  i.e. 

Q(T)  =  a  +  b  T  +  cT2  +dT3 .  (3) 

The  temperature  range  for  the  polynomial  fit  was  determined  by 
requiring  that  the  expression  reproduce  the  calculated  partition  sum 
to  better  than  1%  for  most  of  the  species  in  the  atmosphere.  For 
heavy  systems  with  many  vibrational  degrees  of  freedom  this  is  not 
always  possible,  see  discussion  below.  From  test  calculations  the 
temperature  ranges  in  Table  1  were  adopted. 

The  procedure  is  as  follows.  First  the  total  internal  partition 
sums  are  calculated.  This  is  done  by  one  of  several  methods 
depending  on  the  species  in  question.  If  all  the  energy  levels  are 
known  to  high  values  such  that  the  calculation  of  the  partition  sum 
has  converged,  a  direct  sum  over  the  energy  levels  with  appropriate 
degeneracy  factors  is  made.  When  the  rotational  energies  of  the 


4 


ground  vibrational  state  are  known  but  not  those  of  the  excited 
vibrational  states,  the  product  approximation  is  used.  And  last, 
conditions  for  which  the  energies  are  not  known  to  high  enough 
states  for  the  partition  sums  to  converge,  the  best  available 
analytical  formula  is  used.  The  partition  sums  are  calculated  at  the 
temperature  spacing  given  in  Table  1  and  these  data  are  then  fit  to 
the  polynomial  expression  of  Eq.  (3)  where  the  coefficients  are 
determined  by  fitting  to  the  data  using  a  Simplex  nonlinear 
minimization  algorithm.12  The  fitting  is  based  on  minimizing  the 
sum  of  the  eighth  powers  of  the  percentage  differences  between  the 
model  and  the  calculated  partition  functions.  The  use  of  the 
percentage  differences  reflects  the  way  partition  functions  are  used 
in  calculations,  as  well  as  the  accuracy  of  the  calculations,  better  than 
absolute  differences  would.  The  use  of  least-eighth  powers  as  a 
fitting  criterion  yields  almost  a  minimax  fit.  A  true  minimax  fit 
would  produce  insignificantly  smaller  maximum  differences,  but 
larger  average  differences  in  many  cases.  The  resulting  coefficients 
are  then  stored  in  a  BLOCKDATA  routine  for  recalculating  the  total 
internal  partition  sum  by  a  FORTRAN  program,  TIPS95.FOR,  which  is 
distributed  with  the  HITRAN  database.2 

Several  tests  were  performed  to  insure  the  quality  of  the 
partition  functions.  The  first  test  was  to  calculate  the  TIPS  across  the 
fitting  boundaries  and  at  296  and  300K  to  compare  with  literature 
values.  Calculations  were  made  at  296,  300,  499.99,  500.01, 
1499.99,  1500.01,  2000,  and  2500K.  In  Table  2  the  boundaries  that 
have  error  greater  than  1%  across  the  boundary  are  listed,  these  only 
occurred  at  the  1500K  boundary. 
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Table  1.  Ranges,  temperatures,  grid  spacing  of  fits  and  final  range 
of  application  of  the  polynomial. 


Range 

Temperatures 
used  in  fit 

Temperature 

Grid 

range  of 

polynomial 

1 

70-565  K 

5  K 

70  <  T  <  500 

2 

465-1565  K 

20  K 

500  <  T  <  1500 

3 

1465-3005  K 

20  K 

1500  <  T  <3005 

Table  2  Percent  Error  across  the  Temperature  Range  Boundaries. 


Molecule 

Isotope 

1500K 

boundary 

CD 

38 

1.3% 

02 

66 

1.1% 

02 

68 

2.0% 

o2 

67 

1.4% 

The  calculated  partition  sums  are  also  compared  with  values  found  in 
the  literature.  This  is  limited  to  molecules  and  isotopic  species  that 
have  been  published.  Note  also  that  in  the  comparisons  we  often 
need  to  remove  the  state  independent  nuclear  factor  to  compare  with 
literature  values.  The  reason  for  this  is  that  in  most  applications  the 
ratio  of  the  partition  sums  are  used  thus  the  state  independent 
nuclear  factors  cancel  so  many  authors  leave  this  out.  The  results 
are  presented  in  Table  3  and  notes  on  these  results  are  given  below. 
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h2o- 


co2- 
03  - 

N20- 

co- 

ch4- 

02- 

NO- 

S02- 

N02- 

nh3  - 

HN03  - 


OH  - 
HF  - 


In  addition,  for  181,  171,  and  162,  the  state 

independent  nuclear  factors  are  1,  6,  and  6  and  the  TIPS 
are  179.7,  1074.6,  and  882.95,  respectively.  The  scaled 
values  for  171  and  162  are  179.1  and  147.2 
respectively. 

State  independent  nuclear  factor  removed. 

Possibly  a  sum  over  all  levels  which  is  incorrect  for 
I6O3,  a  Bose  Particle;  Energies  to  J=69  only, 
no  notes, 
no  notes, 
no  notes, 
no  notes. 

neglects  factor  of  3  from  F  states 
I  am  not  sure  why  this  value  is  low. 
no  notes. 

1  am  not  sure  why  this  value  does  not  agree  with  other 
references. 

We  give  the  total  partition  sum  compared  with  Ref.  13 
who  report  only  Qr.  For  most  molecules  there  is  only  a 
small  difference;  however,  for  HNO3  there  are  9 
vibrational  fundamentals  some  of  which  are  low  lying. 
Thus,  the  vibrational  part  contributes  significantly  at 
300K.  The  JPL  value  agrees  well  with  the  classical 
rotational  partition  sum  of  27346.5. 

2  F  levels  per  J  level  omitted 

State  independent  nuclear  factor  removed. 
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HC1  - 


HBr  - 

HI  - 
C10- 


OCS- 


H2CO- 
HOC1  - 
N2  - 
HCN  - 

CH3CI  - 
H202  - 

C2H2- 

c2h6- 


State  independent  nuclear  factor  removed.  We  note 
that  the  classical  partition  sums  are  19.7  for  both 
isotopomers. 

We  note  that  the  classical  partition  sums  are  24.6  for 
both  isotopomers. 

State  independent  nuclear  factor  removed. 

Burkholder  et  al.'s  value14  was  corrected  by  the 
vibrational  partition  sum  (Qv=1.017)  and  multiplied  by 
4  to  account  for  the  F  states.  The  JPL  (Ref.  13)  values  do 
not  have  Qv  contained  in  them. 

Qv  =  1.208  at  300K;  When  this  is  included  with  the  JPL 

value  the  result  is  1241.82,  in  agreement  with  our 

value. 

no  notes. 

no  notes. 

no  notes. 

Values  in  the  table  have  been  multiplied  by  Qv=1.069  to 
compare  with  the  TIPS  value  directly, 
no  notes. 

The  JPL  value  was  multiplied  by  Qv(300K)=l .02  to 
compare  with  the  TIPS  value. 

The  value  of  McDowell  was  multiplied  by 

Qv(300K)=1.186  to  compare  with  the  TIPS  value. 

Both  values  contain  the  state  independent  nuclear  factor 
(gj=64),  the  TIPS  value  is  for  T=296K. 
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ph3  - 


The  JPL  value,  405.64,  has  been  multiplied  by  the 
symmetry  factor  of  McDowell,^  8/3,  and  the  state 
independent  nuclear  factor  (gj=2). 

COF2-  The  value  from  Ref.  23,  Qr=59715,  was  multiplied  by 
QV(296K)  =  .163  to  give  the  value  in  the  table. 

SF6  -  The  TIPS  value  is  at  296K. 

H2S  -  no  notes. 

HCOOH  -  The  JPL  value,  8883.82,  has  been  multiplied  by 

Qv(300)=l .123  and  the  state  independent  nuclear  factor 

(gj=4). 

HO2  -  The  state  independent  nuclear  factor  is  contained  in 

both  values. 

CIONO2-  The  JPL  values  for  the  35C1  and  37C1  species,  100540  and 
103108,  were  multiplied  by  Qv(300)=l  .83 1  and  the 
state  independent  nuclear  factor  ( g j  =  1  2 ) . 

NO+ -  The  state  independent  nuclear  factor  is  contained  in 

both  values. 

Table  3.  Comparison  of  the  TIPS  values  with  literature  values 
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Table  3.  continued 


Temp  K 

Ref. 

Q(Ref) 

O(TIPS) 

12cl6o2 

296 

1  6 

286.11 

286.22 

500 

1  6 

625.76 

625.36 

800 

1  6 

1625.3 

1608.00 

13cl602 

296 

1  6 

288.34 

288.47a 

500 

1  6 

637.93 

638.37 

800 

1  6 

1665.3 

1664.45 

16012C1H0 

300 

1  7 

618.41 

618.47 

16o3 

3  00 

1  3 

6084.16b 

3562.7 

225 

1  3 

4563.52b 

2219. 

296 

18b 

3473 

3481.9 

160180160 

296 

1  8 

3599 

3  645 

160160180 

296 

1  8 

7385 

7462 

14n  ,4N160 

1 

200 

IB 

342.00 

342.9 

300 

mm 

564.08 

564.8 

l2cl6o 

B1 

225 

81.708 

81.72 

300 

mm 

108.787 

108.878 

1  0 


Table  3.  continued 


Temp  K 

Ref. 

Q(Ref) 

O(TIPS) 

12C1H4 

300 

20 

598.639 

602.52 

13C'H4 

300 

20 

598.600 

602.50 

12ClH3D 

300 

20 

1613.818 

1624.27 

16o2 

300 

1  3 

218.655 

218.65 

294 

2  1 

213.9 

214.26 

14N16Q 

3  00 

1  3 

1 159.46 

1 151.38 

296 

22 

380.03b 

1  132 

32sl6o2 

300 

1  3 

5899. 297b 

6497.5 

296 

23 

6345 

6244.5 

14N16o2 

296 

24 

13618 

1  363  1 

14NlH3 

I 

296 

25 

1057. 9b 

300 

26 

580.616 

300 

27 

589.26 

I 

H14Nl603 

300 

1  3 

27938 

35779 

1601H 

220 

28 

27. 7b 

5  8.6 

3  00 

1  3 

81.494 

81.5 

300 

29 

40.75 

81.5 

iH!9f 

300 

30 

39.8 

32.85a 

1H35C1 

300 

1  3 

81.232 

20.35a 

Table  3.  continued 


Temp  K 

Ref. 

Q(Ref) 

Q(TIPS) 

1H37C1 

300 

1  3 

81.352 

20.38a 

IH79Br 

300 

1  3 

101.2045 

25.36a 

1H8,Br 

300 

1  3 

101.2511 

25.36a 

1H127i 

300 

30 

32.025 

32.85a 

35C1160 

296 

■ 

3291b 

3288 

300 

n 

3291.54 

3347.4 

37Q160 

296 

n 

3348b 

3345 

300 

n 

3348.1  1 

3406.0 

16012C32S 

300 

1  3 

1028.25b 

1242 

1H212C160 

300 

1  3 

2876.7 

2742 

1H213C160 

300 

1  3 

2949.7 

2811 

1H16o35C1 

300 

1  3 

2380 

2468 

1H16035C1 

300 

1  3 

2422.5 

2512 

14N2 

298 

3  1 

469.7 

470.28 

1H12C14N 

300 

1  3 

424.154 

mm 

3  00 

6 

151.2b 

fiH 

300 

30 

150.73b 

151.3 

Table  3.  continued 


Temp  K 

Ref. 

Q(Ref) 

O(TIPS) 

12C1H335C1 

300 

30 

14770 

1483  1 

lH2l602 

300 

1  3 

7804b 

7850 

12c2]h6 

300 

8 

55940 

54627b 

3lplH3 

300 

1  3 

3226b 

3322 

12c16o19f2 

296 

32 

69448 

69763 

32S19F6 

300 

7 

1.719  E6 

1.62  E6C 

1H232S 

296 

33 

505.6 

503.2 

1H233s 

296 

33 

506.2 

504.5 

1H234s 

296 

33 

506.9 

503.9 

lH12C160160lH 

300 

1  3 

39906b 

39960 

lH1602 

300 

30 

2191 

2196 

35C116014N1602 

300 

1  3 

2209064b 

2208960 

37cil6o'4N1602 

3  00 

1  3 

2265489b 

2265182 

14n16o+ 

300 

30 

31 1.8 

313.0 

a  State  independent  nuclear  factor  removed, 
b  See  text  for  explanation. 
c  The  TIPS  value  is  at  296K.. 
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We  have  found  that  species  with  many  low  lying  vibrational 
modes  are  not  well  represented  by  the  polynomial  expression.  The 
reason  is  the  vibrational  partition  sum  increases  very  rapidly  with 
temperature  and  thus  the  product  of  the  rotational  and  vibrational 

partition  sums  increases  more  rapidly.  The  4  term  polynomial  has 
difficulty  fitting  the  data  to  the  1%  uncertainty  threshold.  A  number 
of  the  species  do  meet  the  1%  criterion;  all  species  of  CH4  in 
temperature  ranges  2  and  3  and  all  species  of  HNO3  in  all 
temperature  ranges.  This  can  be  understood  by  looking  at  Figures  1, 
2,  and  3  which  show  Q(T)  vs.  Temp  for  H2O,  CH4,  and  HNO3  in  the 
three  temperature  ranges.  HNO3  has  9  vibrational  modes  starting  at 
458  cm*1.  H2O  has  3  vibrational  modes  starting  at  1594  cm-1  and 

CH4  has  9  vibrational  modes  starting  at  1311  cm’1.  In  the  range  70- 
500K  (Figure  1),  H2O  and  CH4  both  rise  slowly  however  HNO3, 
because  of  9  low-lying  vibrational  modes,  is  increasing  rapidly.  H2O 
has  a  fit  error  of  0.37%,  CH4  a  fit  error  of  0.92%,  and  HNO3  a  fit  error 
of  3.6%.  In  Figure  2  the  second  temperature  range  is  shown,  HNO3  is 
going  off  the  scale  and  the  fit  error  is  24%,  H2O  is  still  relatively  flat 

and  the  fit  error  is  0.14%,  however  CH4  is  starting  to  increase  rapidly 

and  the  resulting  fit  error  is  4.9%.  In  the  high  temperature  range  all 
the  vibrational  modes  contribute  to  Q(T)  and  in  Figure  3  we  see  that 
HNO3  is  off  scale  and  CH4  is  going  off  scale  (6.9%  fit  error)  but  H2O  is 
still  reasonable  to  fit  giving  an  error  of  0.03%.  We  must  keep  track 
of  all  species  that  do  not  meet  the  1%  fit  criterion  and  make  a 
decision  on  how  these  will  be  handled  in  the  block  data  BD-QT.FOR 
routine.  At  present,  my  feeling  is  to  flag  these  with  negative 


1500  1750  2000  2250  2500  2750  3000 

Temp 


Figure  3.  Q(T)  for  H2O,  CH4,  and  HNO3  vs.  T  in  the  range  1500- 
3000  K. 

By  studying  the  graphs  of  the  error,  we  determined  that  some  of  the 
coefficients  could  be  used  for  species  where  the  error  was  greater 
than  1%.  Table  4  list  the  species  for  which  the  coefficients  were 
defaulted  (-1,  0.,  0.,  0.)  and  Table  5  list  those  for  which  they  were 
retained  even  though  the  error  was  larger  than  the  criterion. 

The  routine  also  uses  Q(296K).  For  this  the  values  of  the 
partition  sum  at  296  K  were  calculated  from  the  polynomial.  For 
species  where  the  polynomial  fit  was  not  good  (SF6  and  CIONO2)  the 
Q(296K)  values  were  obtained  by  using  an  interpolating  polynomial 
in  Mathematica  with  the  values  of  Q  at  280  to  305  K  in  5K  steps.  The 
values  are  in  Table  6. 
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Table  4.  Species  for  which  Coefficients  are  Set  to  Default  Values 


Species 

ISO 

T  Range 

CH4 

211 

2  and  3 

311 

2  and  3 

211 

2  and  3 

HNO3 

146 

1,  2,  and  3 

CH3CI 

215 

2  and  3 

217 

2  and  3 

C2H6 

1221 

1,  2,  and  3 

COF2 

269 

2  and  3 

sf6 

29 

1,  2,  and  3 

hcooii 

126 

2  and  3 

CIONO2 

5  646 

1,  2,  and  3 

7646 

1,  2,  and  3 

Table  5.  Species  for  which  Coefficients  are  Reported  and  Error>l% 


Species 

ISO 

1 

IFit  Error  1 

2 

3 

H200 

126 

2.2 

2.2 

136 

2.2 

2.2 

128 

2.2 

2.2 

CH3CI 

215 

1.3 

217 

1.3 

PH3 

2.7 

2.5 

COF2 

269 

1.7 

HCOOH 

1  26 

1.7 

Table  6  Q(296K)  from  Interpolating  Polynomial 


§ 

ISO 

0(296K) 

sf6 

29 

1.62242  x  106 

C10N02 

5646 

7646 

2.12829  x  106 
2.18246  x  106 

2.0  SPECTROSCOPIC  PARAMETERS  FOR  THE  OXYGEN  MOLECULE 

Atmospheric  spectra  of  oxygen  are  used  for  deducing 
information  about  properties  of  the  atmosphere  and  other  species  in 
the  atmosphere,  for  example,  measurements  of  O2  emission  are  used 
as  a  standard  for  O3  profile  analysis. 34,35  There  is  a  need  to  have 
available  the  most  accurate  parameters  for  this  molecule.  The 
spectrum  of  the  oxygen  molecule,  even  though  O2  is  a  simple 
diatomic,  is  unexpectedly  complex.  The  oxygen  molecule  has  two 
unpaired  electrons  with  a  total  spin  of  1  in  the  electronic  ground 
state.  There  are  two  low-lying  excited  electronic  states  which  give 
rise  to  near-IR  and  visible  spectra,  and  the  symmetries  of  the 
isotopomers,  16C>2,  160180,  and  160170,  affect  the  number  of  allowed 
states.  In  addition,  both  magnetic  dipole  and  electric  quadrupole 
transitions  occur,  the  degeneracy  of  states  needs  to  be  carefully 
considered,  and  transitions  occur  from  the  microwave  to  the  visible. 
Although  the  line  intensities  are  usually  small,  the  high  mixing  ratio 
and  long  optical  path  in  the  terrestrial  atmosphere  compensate  to 
produce  meaningful  absorption. 3 6 


In  this  work,  the  spectral  parameters  for  the  oxygen  molecule 
are  calculated  for  the  electronic-vibrational  bands  listed  in  Table  7. 
These  data  represent  an  improvement  to  the  data  contained  on  the 
1992  version  of  the  HITRAN  molecular  absorption  database,3  which 
are  from  calculations  made  in  1982.3^  The  calculations  consider  the 
lower  state  energy,  the  wavenumber  of  the  transition,  the  line 
intensity,  and  transition  moment  squared  of  the  spectral  lines.  In 
addition,  halfwidths  as  a  function  of  transition  quantum  number  are 
determined  from  the  available  experimental  measurements.  All 
fundamental  physical  constants  used  in  the  calculations  are  those 
reported  by  Cohen  and  Taylor.38  The  calculations  were  made  in 
double  precision  on  several  different  computer  systems  using  codes 
written  in  FORTRAN.  Below,  we  discuss  the  theory  of  molecular 
oxygen  and  describe  the  improvements  to  the  data  for  each  band. 


2.1  General  Theory 

To  understand  the  spectrum  of  the  molecular  oxygen  in  the 
atmosphere,  one  must  consider  the  properties  of  the  three  most 
abundant  isotopic  species  of  oxygen  and  the  structure  of  the  three 
lowest-lying  electronic  states.  An  excellent  review  for  this  material 
can  be  found  in  Herzberg.5'39 


Table  7. 


Electronic-vibrational  bands  of  O2  in  the  1996  HITRAN 
data. 


Isotopomer 

Electronic  Band 

1602 

160180 

160170 

V’  4r- 

v” 

v’  4-  v” 

v’  4-  v” 

X3z-g  4-  X3l-g 

0  <- 

0 

0  4-  0 

0  4-  0 

1  4- 

0 

1  4- 

1 

ti^Ag  4 — 

0  4- 

1 

0^0 

0  4- 

0 

1  4- 

0 

b^L^g  4 —  X3Z'g 

0  <- 

0 

0  4-  0 

1  4- 

0 

1  4-  0 

1  4-  0 

2  4- 

0 

2  4-0 

1  4- 

1 

0  <- 

1 

The  theoretical  treatment  of  the  ground  state  of  molecular 
oxygen  was  first  given  by  Tinkham  and  Strandberg411  and  it  was  later 
clarified  with  the  aid  of  the  transformation  theory  of  spherical 
tensors  by  Steinbach  and  Gordy.41>42  References  43  and  44  give 
excellent  overviews  of  the  theory.  The  ground  electronic  state  of  the 
oxygen  molecule  is  the  X^Z'g  state.  There  are  also  two  low-lying 
electronic  states,  the  a^Ag  state  and  the  b'Z+g  state  at  7918  cm4  and 
13195  cm-1  above  the  ground  state,  respectively.  Since  the 
electronic  ground  state  is  the  X3£-g  state,  the  average  electronic 
orbital  angular  momentum  vanishes  (A=0).  Still,  there  is  an 
instantaneous  non-zero  value  of  the  orbital  angular  momentum 
which  produces  a  precessing  magnetic  dipole  moment  of  orbital 
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origin.  The  total  electronic  spin  is  S=l,  so  that  molecular  oxygen  has 
a  permanent  magnetic  dipole  moment  of  approximately  two  Bohr 
magnetons.  The  energy  is  described  in  terms  of  the  total  electronic 
spin  vector,  S,  and  the  rotational  angular  momentum,  N.  In  the 
electronic  ground  state,  Hund's  coupling  case  (b)  dominates39  (more 
recent  work  shows  the  need  for  intermediate  coupling)  and  the  total 
angular  momentum  J  is  J=N+S.  Thus,  for  each  value  of  N  there  are  3 
allowed  values  of  J.  The  oxygen  molecule  can  interact  with  an 
electromagnetic  field  via  its  magnetic  dipole  (md)  moment  or  by  its 
electric  quadrupole  (eq)  moment.  In  this  work,  calculations  are 
made  for  both  md  and  eq  lines  of  O2. 

The  principal  isotopic  species,  16C>2,  is  unique  in  that  it  is  the 
only  homonuclear  diatomic  molecule  of  the  three  isotopomers  on  the 
database.  Because  of  this  symmetry  and  the  fact  that  the  oxygen-16 
nuclei  have  zero  nuclear  spin,  the  molecule  behaves  as  a  Bose 
particle.  Thus,  the  total  wavefunction  must  be  symmetric  with 
respect  to  inversion  through  the  center.  The  spins  of  the  individual 
electrons  form  a  resultant  S=1  that  gives  rise  to  a  coupling  of 
electronic  (S)  and  rotational  angular  momentum  (N)  yielding  a  triplet 
of  states  labeled  by  J;  J  =  N,  J  =  N  +  1,  and  J  =  N  -  1.  The 
wavefunctions  for  the  vibrational  and  nuclear  motion  are  symmetric, 
whereas  the  electronic  wavefunction  for  a  L'g  state  is  antisymmetric, 
indicating  the  only  allowed  rotational  wave  functions  are  those  for 
odd  (antisymmetric)  N.  The  result  is  that  one-half  (all  N  even)  of  the 
states  of  1602  are  missing  in  the  ground  electronic  state. 

In  the  a1  Ag  electronic  state,  the  resultant  of  the  individual 
electron  spin  is  zero  (S=0)  giving  J=N,  but  the  rotational  levels  are 


split  into  two,  one  symmetric  (+)  and  one  antisymmetric  (-)  state,  by 
A -doubling.  Because  of  the  constraints  placed  on  the  wavefunction 
by  Bose-Einstein  statistics,  only  the  (+)  state  is  allowed  and,  although 
all  N  are  occupied,  only  one-half  of  the  states  are  realized.  The  b^+g 
electronic  state  also  has  S=0  giving  N=J  and  the  electronic 
wavefunction  for  the  I+g  state  is  symmetric  as  are  'Fnuc  and  'Fyib; 
thus,  the  only  rotational  wavefunctions  are  the  symmetric  (N  even) 
ones.  Because  of  the  symmetry  of  the  1 6C>2  species  and  the  nuclear 
spin  I(16O)=0,  one-half  of  the  rotational  states  of  each  electronic  level 
are  missing. 

For  the  two  isotopic  species  on  the  database,  16()18o  and  ^O^O, 
inversion  through  the  center  is  no  longer  a  valid  symmetry 
operation.  For  this  condition,  all  rotational  levels  of  the  molecule  are 
allowed,  both  (+)  and  (-).  The  ground  electronic  state  is  a  triplet,  J=N 
and  J=N±  1,  with  all  N  allowed,  N=l,2,3,«**.  The  two  low-lying 
electronic  states,  a^Ag  and  b^Z+g  are  singlets  (J=N)  with  all  values  of 
N  (even  and  odd)  starting  at  Nmjn=  A.  Note  that  for  the  a^Ag  state 
there  is  A -type  doubling,  thus  each  N  has  two  states. 

Figure  4  is  a  rotational  state  diagram  for  the  lowest  three 
electronic  states  of  O2.  A  note  of  caution  here:  the  rotational  states  of 
a  diatomic  molecule  are  classified  according  to  the  behavior  of  the 
total  wavefunction  with  respect  to  reflection  at  the  origin  and  not  of 
the  rotational  wavefunction  alone  (see  pp  128-129  of  Ref.  39).  Thus 
in  the  diagram,  the  states  labeled  (+)  in  the  X^Z'g  state  correspond  to 
N=odd  (antisymmetric)  states;  whereas  the  states  labeled  (+)  in  the 
b^+g  state  correspond  to  N=even  (symmetric)  states.  The  principal 
species  of  oxygen,  ^02>  follows  Bose-Einstein  statistics  and  as  such 
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does  not  have  any  of  the  (-)  symmetry  levels.  Thus,  in  the  diagram 
all  (-)  states  are  eliminated,  in  the  X3Z'g  state  all  even  N  levels  are 
missing,  in  the  a]Ag  state  one  level  for  each  N  is  missing,  and  in  the 
b1£+g  state  all  odd  N  are  missing.  For  the  non-homonuclear  species 
all  levels  are  allowed.  These  symmetry  factors  have  an  important 
role  in  the  statistics  of  the  molecule.  Taking  the  ratio  of  the  number 
of  levels  in  the  X3Z-g  state  to  the  a!Ag  state  gives  3:2  regardless  of 
the  isotopomer  (this  is  because  for  16C>2  one-half  of  the  levels  are 
missing  in  each  electronic  state).  Likewise  the  similar  ratio  of  the 
X3I*g  state  to  the  b1L+g  state  is  3:1. 
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Figure  4.  State  diagram  for  the  3  lowest  electronic  levels  of  the 

oxygen  molecule,  X3Z’g)  aUg,  and  b^+g,  and  their 
statistical  weights. 

2.2  Line  Intensity 

Various  formulas  are  used  to  calculate  the  line  intensities 
depending  on  the  particular  transition.  Below  these  formulas  are 
presented  along  with  other  relationships  needed  for  the  calculations. 
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Also  presented  are  relationships  to  explain  problems  encountered  in 
determining  the  Einstein-A  coefficients  for  transitions  in  the  X^Z'g  — » 
a1  Ag  band.  The  transitions  involve,  in  general,  different  electronic, 
lambda  doubling,  vibrational,  and  rotational  states  which  will  be 
labeled  by  e,  A,  v,  and  J,M,  respectively  with  the  usual  spectroscopic 
notation  of  a  single  prime  for  the  upper  state  and  a  double  prime  for 
the  lower  state.  In  the  absence  of  an  electric  or  magnetic  field,  the  M 
states  are  degenerate  and  will  be  summed  over  to  express  the  line 
intensity  for  the  transition  labeled  by  e"A  "v"J"— >e'A  VJ’.  In  the 
following,  a  shorthand  notation  will  be  adopted  by  substituting  T]  for 
the  set  of  quantum  numbers  eAvJ.  Thus,  q"  =  e"A"v"J"  and  t|'  =  e'A'v'J’. 
The  line  intensity  in  units  of  (cm- l/molec*cm‘2)  for  a  transition  tj"  — > 
rj '  is  (see  Refs.  1  and  45) 

c  —  _8jlL  i  (<V  e'E^~/kT)  VqvxtV) 

3 he  a  dr  Q(T)  c 

x  (l-  e'hvcn"K,iVkT)  |R(£"a,v'J"M")(£'a'vTM'J2  10  36,  (4) 

M",  M' 


where  the  total  internal  partition  sum,  Q(T),  is 

Q  (T)  =  X  de  dA  dv  dj  dsym  e-E^v j/kT  ,  (5) 

EAvJ 

and  Ia  is  the  isotopic  abundance  of  the  species,  the  factor  10*3 6  is 
needed  for  the  units  chosen,  IR  (e"A"v"J"M")(e’A  VJ'M')I2  is  the 
transition-moment  squared  and  has  units  of  Debye  squared/molec, 
d^"  is  the  degeneracy  and  E^"  the  energy  of  the  state  r|",  V(t|")(ti’)/c  = 
(o (-p " >(t| ')  is  the  wavenumber  of  the  transition,  and  all  other  symbols 
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are  the  usual  constants.  (Note,  in  Eq.  (4)  and  in  a  few  expressions 
that  follow  we  have  not  canceled  some  of  the  degeneracy  factors  in 
order  to  emphasize  the  fractional  number  of  molecules  in  the  state 
Nij ".)  In  this  report  we  use  the  "regular:  matrix  elements  while  in 
the  Gamache  and  Rothman  1992  paper1  we  used  the  "weighted"  ones. 
The  degeneracy  factors  given  by  the  product  of  the  degeneracy 
factors  for  the  quantized  motions,  dT1=dedAdvdjdSym  with  de=(2S+l), 
dA=(2-5A,o),  dv=l,  dj=(2J+l ),  and  dsym  is  one  for  the  heteronuclear 

species  and  the  (+)  states  of  the  homonuclear  species  and  zero  for  the 
(-)  states  of  the  homonuclear  species.  This  factor  accounts  for  the 
symmetry  restrictions  for  the  homonuclear  diatomics  (note,  this 
gives  an  overall  factor  of  1/2  in  the  partition  sum). 

This  expression  can  be  written  in  terms  of  the  Einstein-A 
coefficient  by  making  use  of  the  relationship4 


_  64k4  IQ’36  ^3  J_  V 
?  ,  w(n  )bi )  , 

3  h  dj'  M-\M' 


|R(E"A"v"J"M")(e•AVJ•M,: 


whence  the  line  intensity  is 


S 


=  la ^ 11 — 1 - 1  — 1 - (l-  e*hV(’1")(,1':,/,kT)  . 

8ttc  Q(T)  dj" 


(7) 


The  Einstein-A  coefficient  is  usually  determined  from  a 
measurement  of  a  vibrational  band  intensity,  Sg.^.y.  so  it  is  useful  to 
formulate  Eqs.  (4)  and  (7)  in  terms  of  this  quantity.  The  integrated 
intensity  for  an  electronic-vibration-rotation  band  in  units  of  cm 
Vfmolec'cm-2)  is  given  by45 
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„eMAMv 

^e'AV 


l(e"A"v")(E,AV) 


R(E"A"v")(e’AV)|  •  (8) 


_  8n3  10'36  m 
3hc 


(dE-A’y  e~Ec"A"v'/kT) 
QeAv 


(Note,  this  is  similar  to  <x/Nl  (Ref.  45,  pi  53  with 
a  =  X  Se«a vr-KAVj' )  with  the  exception  that  the  radiation  field 
J"J' 

term,  (i-  e-hv(t)")(Ti')/kT^  explicitly  retained  in  the  rotational  part  of  the 
intensity  formula.  Thus,  our  definition  of  does  not  contain  the 

approximate  factor  (l-  e‘hv<c'A"v''><e'AV>/kT).)  To  implement  Eq.  (8),  the 
electronic-vibration  and  rotation  parts  are  separated  using  the 
product  approximation  for  the  total  partition  sum  and  the  transition- 
moment  squared,  and  the  assumption  of  additivity  of  energy 


Qtot  —  QeAv*Qrol  > 

|R(t|")(t|')|  =  |R(e"A'V')(e'AV)|  •|R(J")(J'J  ’ 

M”,M’ 

Et|"  =  Eema"v"  +  Ejm  . 


Inserting  these  into  Eq.  (4)  and  rearranging  terms  yields 


S,.^.  -  ^  '»-3<  C'l  ‘  '  '  ,.rA'  ,'R.r-AV.fAvf 

3  h  c  QeAv 

X  Ia  (dr  dSym  e'Ej  /k'[)  — —  ( 1 .  e-hv<T'">w>/kT)  |R(d(J'J2  .  (10) 

dj"  djjyjn  Qrot  “(e  ^  v  )(cAv) 


The  first  part  of  this  expression  is  simply  the  band  intensity  as 
defined  in  Eq.  (8);  thus  the  line  intensity  can  be  written 
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S-n" 


n  -»ri 


=  4Av”  Ia  (dj"  dsym  e^J,'/kT)  _ (i_  e'hv(T>"><ii')/kT) |R(j-)(jJ2 


dj"  dsynl  Qrot 


0)(£"A"v")(£'AV) 


(11) 


In  what  follows  we  will  need  the  relationships  between  the  Einstein 
coefficients  and  the  electronic-vibrational  transition-moment 
squared.  They  are 


R . .  8 7t3  10'36  Ip  ,.A„  „v  ...  J2 

*J£  A  v  — >eAv  —  |K‘(£  A  v  )(eAv]  i 


3h‘ 


Brv.--.  A-.-  =  87,3  |R<e-a  v-xvavii  , 

3h‘ 


dEAV 


(12a) 

(12b) 


and 


-  64k*  10'36  d 


A,  „  „  _  W  "T  ^  t  1  U _ 

e’A  v  — >e  A  v  —  uj/n’Av 

3  h 


(e'A'v')(e"A"v“)  ,E  A  V  |R(e"A”v")(£'AvJ  •  (12c) 

de’AV 


Solving  Eq.  (12c)  for  the  transition-moment  squared  and  inserting 
into  Eq.  (8)  gives  the  result 


<j£"A"v" 

^e'A'v' 


1  (dE"A”v"  e~E£"A"v"/kT) 

QeAv  W(£'AV)(£''A"v") 


de'A’v'  ^ 
de"A'V' 


e'A 


V— >e"A"v 


(13) 


This  expression  allows  the  Einstein-A  coefficient  to  be  determined  by 
measuring  the  integrated  band  intensity.  This  relationship  with  Eq. 
(11)  allows  the  e"A "v"J"->e'A 'v'J'  line  intensity  to  be  written  in  terms 

of  the  Einstein  Ae'AV->e"A"v"  coefficient 


S 


la 

&KC 


(d 


£"AVT 


dr  Qtot 


^(E,,A,V,J,,)(£,A,v,J,) 

^(e"A"v”)(e'AV) 


27 


(14) 


X  (l-  e-hvcnw^T)  Lj.  Ae.AV^e..A.v. 

de"A"v" 

where  the  partition  function  and  the  energy  are  no  longer 

approximated  and  the  Honl-London  factor  is  inserted  for  the 
rotational  transition-moment  squared.  Note  that  the  degeneracy 
factors  preceding  the  Einstein-A  coefficient  account  for  the  number 

of  J  and  M  states  in  the  upper  and  lower  electronic-vibration  states 
and  the  degeneracy  factors  with  the  partition  function  term  are 
deAvJ=(2S  +  l)(2-5A>o)(dSyin)(2J  +  l).  From  the  symmetry  arguments 
presented  above,  it  is  clear  that  the  ratio  of  the  degeneracy  factors, 

de’A v/de'-A'W,  is  2/3  for  X3Z'g  -»  aJAg  transitions  and  1/3  for  the 

X3X'g  -»  b1Z+g  transitions. 

For  several  of  the  bands  calculated,  the  program  employs  the 
Einstein-A  coefficients  in  units  of  sec‘1.  Often  one  must  work  from 
measured  vibrational  band  intensities.  In  order  to  make  use  of  these 
values,  we  must  convert  from  vibrational  band  intensity  to  the 
Einstein-A  coefficient.  This  relationship  was  presented  in  Eq.  (13). 
However,  quite  often  in  the  literature  the  vibrational  band  intensity 
is  reported  in  units  of  crtr'  kirr^  atm‘1  STP.  The  values  must  be 
converted  to  the  HITRAN  units  of  cm'1/(molec*cm-2)  to  apply  the 
equations  presented  here.  This  is  accomplished  by  the  relationship 

sM _ cm~' - 1 x — Nk — x  l05£JQ-=Sv'(cm‘1  km'1  atm'1  STP)  (15) 

V  'molcc  cm'2'  p(atm)  km 

where  Nl  is  Loschmidt's  number,  the  number  of  molecules  per  cubic 
centimeter  of  perfect  gas  at  STP.  Thus,  we  have 
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Sy.  ( - smd - )x  2.6867*  1 024-c-m-=  Sy-  (cm'1  km’1  atm1  ST p)  .  (16) 

'molec  cm'2'  k  m 


From  Eq.  (13),  the  Einstein-A  coefficients  used  in  the  program 
have  been  updated  using  vibrational  band  intensities  reported  in  the 
literature.  The  vibrational  band  intensities  and  Einstein-A 
coefficients  used  in  the  program  are  presented  in  Table  8  for  the 
b^+g  <-  X3£-g  band. 


Table  8.  Vibrational  band  intensities  and  Einstein-A  coefficients 
for  the  b^+g  <—  X3S‘g  band  of  O2. 


Band 

mm 

Sv.  b 

°V 

A  c 

Ref. 

A(ref.)c 

A(02CALC)C 

A  band 

532. 

1.9812-22 

0.0770 

IBB 

0.077 

0.077 

582. 

2.1712-22 

0.084 

47 

0.084 

NA 

2.2812-22 

0.0887 

2 1 

0.0887 

NA 

B  band 

40.8 

1.5212-23 

0.00724 

48 

NG 

0.00591  * 

38.8 

1.4412-23 

0.00689 

47 

0.0069 

NA 

1.52 

5.6612-25 

3.2412-4 

eh 

NG 

2.212-4* 

RflSfl 

1.50 

5.5812-25 

3.1912-4 

47 

3.212-4 

NA 

1.26 

4.7112-25 

2.6912-4 

50 

NG 

NA 

0.0269 

1.0012-26 

6.7312-6 

47 

6.712-6 

NA 

RBI 

0.136 

5.0612-26 

0.0701 

47 

0.0704 

0.0704 

mm 

0.0114 

4.2412-27 

0.00467 

47 

0.0047 

0.0047 

a  units  of  cm-1  km-1  atm'1  STP 
b  units  of  cm'1/(molec  cm"2) 
c  units  of  sec’Vmolec 

*  does  not  agree  with  literature  value  (see  text) 
NA  not  applicable 
NG  not  given 
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The  vibrational  band  centers  are  computed  using  the  following 


expressions 

G(v)  =  G>e  (v+1/2)  -  (0exe  (v+1/2)2  +  (Oeyc  (v+1/2)3  -  C0eze  (v+1/2)4  ,(17) 

G(v’)  =  to'e  (v’+l/2)  -  0)ex^  (V+1/2)2  +  0)eye  (V+1/2)3  -  0)eZe  (v'+l/2)4  ,(18) 
and 

CD0  =  Tc  +  G(V)  -  G(v)  .  (19) 

The  constants  (Te,  coe,  coexe>  ioeye,  toezc,  a>c\  coeXe',  toeye’,  <*>eZe')  are  taken 
from  Krupenie51  for  both  the  a1  A+g  and  b^+g  electronic  states. 

The  partition  sums  used  in  the  program  are  those  from  the 
TIPS52  (Total  Internal  Partition  Sum)  program.  Before  addition  to  the 
HITRAN  database,  the  line  intensities  are  filtered  through  a 
wavenumber  dependent  cutoff  given  by 


\1  +  e  kTl 


(20) 


with  o)c=2000  cm*1  and  S(oc=3.7  x  10’30  cnr  Vtmolec^cnr2).  All  lines 
with  an  intensity  less  than  the  cutoff  are  not  included  in  the 
database  with  the  exceptions  noted  below. 
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2.3  Halfwidths 


Measurements  of  the  halfwidths  are  available  for  several 
bands  of  the  O2  molecule.  The  data  of  Krupenie51  are  taken  for  the 
X3I’g  pure  rotation  bands,  the  A-band  data  are  from  Ritter  and 
Wilkerson,21  the  B-band  data  from  Giver  et  al.,^8  and  the  y-band  data 
are  from  Melieres  et  al.50  The  data  of  Refs.  21,49-50  are  plotted  in 
Figure  5  as  a  function  of  N"  for  the  various  types  of  transitions.  The 
transitions  are  labeled  by  ANAJ  N"  J";  however  with  symmetry 
arguments  the  notation  aNaJn"  is  often  used.  For  the  A-band  we  find 
that  if  the  RR  and  RQ  values  are  shifted,  (N"  N”+2),  they  agree  with 

the  pP  and  PQ  values  (see  Figure  6),  i.e. 


y(ppN")  =  y(RRN"+2) 


and 


(21) 


7(pQn")  =  Y(rQn-+2) 


For  the  B-band,  the  difference  between  the  pP,  PQ,  RR,  and  RQ  values 
for  a  given  N"  is  small.  The  data  of  Melieres  et  al.^O  for  the  y-band, 
however,  show  large  variations  in  the  halfwidth  as  a  function  of  N" 
not  seen  in  the  other  data  (see  open  circles  and  triangles  in  Figure  5). 


PP-A 


*» -  PQ-A 


Figure  5.  Halfwidth  data  cnr1  atm1  at  296K  for  the  A-,  B-  and 
y-bands  of  O2  according  to  branches  as  a  function  of  N. 


Comparing  the  values  from  Refs.  21,49-50,  we  find  that  the  A- 
and  y-band  results  agree  with  each  other  and  the  B-band  results  are 
some  10%  larger.  The  data  selected  for  O2  on  the  1996  HITRAN 

database  use  the  following  procedure  for  the  halfwidths.  The  X^L'g 
pure  rotation  band  uses  the  data  reported  by  Krupenie^1  for  the  60 
GHz  lines.  The  electric  quadrupole  transitions  and  the  transitions 

involving  the  a]Ag  state  use  the  A-band  values.21  For  the  A-  and  y- 
bands,  the  halfwidths  of  Ritter  and  Wilkerson21  are  averaged  as  a 
function  of  N"  for  the  pP  and  PQ  lines  together  and  the  RR  and  RQ 
lines  together.  We  have  calculated  transitions  to  N"=80  and  the 
measurements  end  at  N"=29;  thus,  it  is  necessary  to  extrapolate  the 
data.  From  the  plots,  an  asymptotic  limit  of  y=0.032  cm-1  /atm  for 
N">40  is  estimated.  Both  the  R  and  P  lines  use  the  same  asymptotic 
limit.  The  values  of  the  halfwidths  for  the  A-  and  B-bands  of  O2  used 
in  HITRAN96  are  given  in  Table  9. 

2.4  Updates  to  the  Oxygen  Data  for  the  1996  HITRAN  Database 

Below,  the  changes  made  to  the  data  for  the  1996  HITRAN 

database  are  discussed  for  each  band  of  molecular  oxygen 

considered.  The  line  position  and  energy  differences  are  defined  as 
the  HITRAN92  value  minus  the  HITRAN96  value.  For  the  line 
intensities  the  ratio  computed  is  HITRAN92/HITRAN96.  The  average 
and  maximum  differences  were  calculated  for  the  line  position  and 
energy  parameters.  For  the  line  intensities,  the  average,  maximum 
and  minimum  ratios  were  calculated. 
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Table  9.  Halfwidths  in  cnrVatm  for  the  A-  and  B-bands  of  O2. 


N" 

A-Pf 

A-Rft 

bitf 

N" 

A-Pf 

A-Rft 

Blit 

1 

.0592 

.0587 

.0616 

2  1 

.0417 

.0406 

.0454 

2 

.0574 

.0562 

.0600 

22 

.0409 

.0402 

.0447 

3 

.0557 

.0538 

.0584 

23 

.0402 

.0397 

.0440 

4 

.0544 

.0523 

.0571 

BM 

.0395 

.0386 

.0436 

5 

.0531 

.0509 

.0558 

25 

.0389 

.0376 

.0432 

6 

.0520 

.0503 

.0555 

26 

.0380 

.0371 

.0428 

7 

.0509 

.0498 

.0552 

27 

.0371 

.0366 

.0422 

8 

.0502 

.0492 

.0547 

28 

.0365 

.0183 

.0418 

9 

.0495 

.0486 

.0541 

29 

.0358 

.0358 

.0414 

1  0 

.0488 

.0479 

.0531 

30* 

.0353 

.0353 

.0410 

1  1 

.0482 

.0473 

.0521 

31* 

.0349 

.0349 

.0407 

1  2 

.0475 

.0467 

.0515 

32* 

.0345 

.0345 

.0404 

1  3 

.0468 

.0461 

.0509 

33* 

.0340 

.0340 

.0400 

1  4 

.0464 

.0455  : 

.0504 

34* 

.0337 

.0337 

.0398 

1  5 

.0459 

.0449 

.0499 

35* 

.0334 

.0334 

.0395 

1  6 

.0452 

.0441 

.0488 

36* 

.0330 

.0330 

.0394 

1  7 

.0445 

.0432 

.0477 

37* 

.0328 

.0328 

.0393 

1  8 

.0438 

.0426 

.0469 

38* 

.0324 

.0324 

.0391 

1  9 

.0431 

.0421 

.0461 

39* 

.0322 

.0322 

.0390 

20 

.0424 

.0413 

.0457 

40* 

.0320 

.0320 

.0390 

A-band  pP  and  PQ  transitions  *  Extrapolated 


t  A-band  RR  and  RQ  transitions 

Iff  B-band  transitions 
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2.4.1  The  principal  isotopic  species,  16C>2. 

The  X^L-g,  (v=01  ±z  (v=01  band 

The  energy  levels  for  the  vibrational  ground  state  of  the  X3I'g 
electronic  state  of  1602  are  calculated  using  the  formalism  of  Rouille 
et  al.53  In  this  work,  the  Hamiltonian  included  all  rotational  terms  to 
second  order54  and  some  terms  to  third  order.55*56  The  molecular 
constants  are  those  of  Rouille  et  al.53  These  are  compared  with  the 
energy  values  from  the  previous  database57  in  Figure  7.  The 
difference  in  energy  between  the  two  formulations  is  given  as  a 

function  of  N  and  J.  For  the  vibrational  ground  state  the  difference  is 

near  zero  for  N  up  to  35  then  rapidly  goes  to  -0.3  cnr1  at  N=80. 

While  the  transition  frequencies  from  this  formulation  differ  only 
slightly  from  the  previous  results,  0.0039  cnr1  maximum  difference, 
7  x  10'5  cm-1  on  average,  the  formalism  of  Rouille  et  al.53  gives 
better  agreement  with  measurements.5  ^ 

This  band  has  both  magnetic  dipole  and  electric  quadrupole 
transitions.  For  the  electric  quadrupole  transitions,  a  newer  value  of 
the  quadrupole  moment  derived  from  far-IR  PIA  spectra, 5^  0.34  x 
10'26  esu  cm2,  has  been  adopted.  There  are  no  data  to  validate  the 
intensities  we  obtain.  With  this  value  for  the  quadrupole  moment, 

the  line  intensities  and  transition-moments  squared  are  a  factor  of 
5.8  weaker  than  previous  calculations.611  Because  of  this  reduction  in 
the  line  intensities,  no  electric  quadrupole  lines  survived  the  cutoff 
for  the  1996  data  set.  This  is  further  discussed  below  for  the  0<-  1 
vibrational  band  of  this  electronic  band.  The  1996  calculations  of  the 
magnetic  dipole  intensities  are  on  average  4%  stronger  than  the  1992 


values. 


This  is  due  to  improved  energy  formulation  and  partition 


sums. 


Figure  7.  Energy  difference  in  cm'1  between  the  formalism  of 
Rouille  et  al.53  and  Ref.  57  versus  N  for  the  v"=0  band 
of  X3X-g  of  16C>2. 


The  (v=l)  <-  XlZ-g  (v=0)  band 

The  1992  HITRAN  database  contained  only  electric  quadrupole 
(eq)  lines  for  this  band.  The  newer  data  contain  both  magnetic 
dipole  (md)  and  electric  quadrupole  (eq)  transitions.61  Thus,  the 
comparison  made  here  is  for  the  electric  quadrupole  lines.  The 
energy  levels  for  the  vibrational  ground  (v=0)  and  first  fundamental 
(v=l)  of  the  X3£-g  electronic  state  of  16C>2  are  calculated  using  the 
formalism  and  constants  of  Rouille  et  al.63  These  energies  are 
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compared  with  the  values  from  the  previous  database37  in  Figure  8 
for  the  v"=l  band.  The  difference  is  small  for  N  up  to  ~20  then 

quickly  goes  to  roughly  5  cm"1  at  N=80.  For  the  transitions  that 

make  the  cutoff  criterion,  N<31,  the  maximum  difference  in  the 
energy  values  is  0.0024  cm'l,  0.0005  cirri  on  average.  The 

corresponding  average  and  maximum  differences  in  the  line 
positions  are  -0.0104  and  0.0850  cm-1.  The  electric  quadrupole 
transition  line  intensities  are  calculated  using  the  formalism  of 

Goldman  et  al.61  The  absolute  intensities  are  determined  by  scaling 
the  relative  intensities  to  the  measurements  of  Reid  et  al.6  2 
Comparison  of  these  intensities  with  the  1992  HITRAN  values  implies 
a  electric  quadrupole  moment  of  0.145  x  10-26  esu  cm2,  not  in 
agreement  with  the  quadrupole  moment  derived  from  far-IR  PIA 
spectra^  (see  previous  section).  This  is  under  investigation. 
Comparing  the  1996  and  1992  values  for  the  intensities  we  find  an 
average  ratio  of  1.0027,  with  some  ratios  between  0.740  and  1.31. 
These  large  ratios  occur  for  33  out  of  the  146  lines  and  occur  only  for 
forbidden  (weaker)  lines,  with  AN*AJ  and  at  low  J.  . 

The  1996  calculations  for  this  band  are  not  filtered  through  a 
cutoff  procedure  and  many  of  the  lines  will  have  very  small 
intensities.  (In  fact,  one  zero  intensity  magnetic  dipole  and  electric 
quadrupole  transition  have  been  retained  in  the  data  for  theoretical 
considerations;  it  helps  to  see  the  effects  of  assumed  parameters  on 
these  lines.)  There  are  254  md  transitions  and  183  eq  transitions  for 
this  band. 
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Figure  8.  Energy  difference  in  cnr1  between  the  formalism  of 
Rouille  et  al.53  and  Ref.  57  versus  N  for  the  v"=l  band 
of  X3Z-g  of  ,602. 

The  fv=n  <-  XZ-l-g  (v=l)  band 

As  discussed  above,  the  energy  differences  from  the  previous 
calculations  and  the  formulation  of  Rouille  et  al.53  for  the  v=l  level 
approach  5  cm-1  at  N=80.  The  maximum  energy  difference  in  the 
lines  in  the  1996  data  is  0.750  cnr1  and  the  average  difference  is 
-0.0953  cm'1.  The  maximum  difference  in  line  positions  is  0.0457 
cm-1  with  the  average  difference  being  -0.00280  cm'1.  The  intensity 
ratios  are  between  0.994  and  0.988  with  an  average  ratio  of  0.990. 
Most  of  this  difference  comes  from  using  an  improved  partition  sum 
in  the  calculations. 
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The  alAg_  (v=01  <-  X$L-&  (v=01  band 

There  are  wavenumber  differences  which  arise  from  a  change 
in  the  energy  formulation  of  X3£'g  (v=0)  to  that  of  Rouille  et  al.53  and 
for  aJAg  (v=0)  to  that  of  Scalabrin  et  al.63  with  the  constants  of  Hillig 
et  al.64  The  maximum  difference  is  0.0160  cm'1  at  N"=37  and  the 
average  difference  is  -0.00234  cm-1. 

The  line  intensities  for  this  band  are  calculated  using  Eq.  (14). 
There  are  three  measurements  of  the  band  intensity  in  the  literature 
(scaled  to  296  K)  which  differ  by  a  factor  of  about  4.5.  The 

measurement  of  Badger  et  al.65  gives  Soo  =  3.6  x  10' 2  4 
cm'1/(molec*cm'2).  The  measurement  of  Lin  et  al.66  is  Soo  =  9.4  x  10' 
24  cm'1/(molec*cm'2),  and  Hsu  et  al.67  report  a  value  of  Soo  =  2.1  x 
10*24  cm"1/(molec*cm*2).  From  these  values  the  authors  determine 
the  Einstein-A  coefficient.  Unfortunately,  some  of  the  authors  have 
used  incorrect  statistical  degeneracy  factors,  d//du;  Badger  et  al.65 

used  3/2  whereas  Lin  et  al.66  and  Hsu  et  al.67  use  3/1.  It  was 
demonstrated  above  that  in  Eq.  (14)  d/=3  and  du=2  .  Thus, 

application  of  Eq.  (14)  (or  Eq.  (13))  will  only  generate  consistent  line 
intensities  if  the  statistical  degeneracy  factors  and  the  derived 

Einstein-A  coefficient  are  from  the  same  author.  Because  of  the  large 
discrepancy  between  the  measured  band  intensities,  we  sought  to 
determine  the  band  intensity  from  another  source.  We  have 
obtained  unpublished  line  intensity  measurements6 8  of  14 

transitions  in  the  a'Ag  (v=0)  <—  X3Z'g  (v=0)  band  that  were 

mentioned  in  the  work  of  Wallace  and  Livingston. 6^  Using  our  O2 
program,  the  Einstein-A  coefficient  used  to  calculate  the  line 
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intensities  was  scaled  to  match  the  measurements  of  Ref.  68.  The 
results  are  in  Table  10.  This  fit  gives  a  band  intensity  of  Soo=3.69  x 
10*24  cm* l/(molec*cirr2)  which  is  close  to  the  Badger  et  al.65  value. 
The  Einstein-A  coefficients,  the  band  intensities,  and  the  statistical 
degeneracy  factors  as  related  by  Eq.  (13)  are  listed  in  Table  11.  The 
calculation  of  the  line  intensities  for  the  1996  database  used  the 
Einstein-A  coefficient  A=2.59  x  10*4  sec*1  with  statistical  degeneracy 
factors  of  d/=3  and  du=2.  The  resulting  line  intensities  are  larger 

than  HITRAN92  by  roughly  a  factor  of  2,  this  is  due  to  incorrect 
inversion  from  Badger  et  al.'s  A  to  S  used  in  previous  versions  of 
HITRAN.  Suspicion  of  such  missing  factors  of  2,  and  concerns  about 
the  interpretation  of  upper  atmosphere  emissions,  such  as  inferring 
ozone  from  SME  (Solar  Mesosphere  Explorer  via  the  a  *Ag  1.27pm 
airglow)  were  expressed  by  Mlynczak  and  Nesbitt7^  (who 
conjectured  a  significant  change  in  A  from  one  of  the  reported  S 
values  in  Table  5  (Hsu  et  al.))  and  by  Pendelton  et  al.71  Recent 
observations  in  the  mesosphere72  confirm  the  Badger  et  al.65  value  of 
the  Einstein-A  coefficient.  Moreover,  new  line  intensity 
measurements  made  at  National  Institute  for  Standards  and 
Technology73  and  at  Rutherford  Appleton  Laboratory74  are  roughly 
15%  larger  than  the  HITRAN96  values.  Preliminary  comparisons 
indicate  agreement  between  these  two  new  independent  high- 
resolution  studies.  When  completed,  these  data  will  be  incorporated 
into  the  next  edition  of  HITRAN. 
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Table  10. 


Measured  intensities  (Brault  and  Brown68)  of  the  a*Ag 
(v=0)  X3I‘g  (v=0)  band,  final  HITRAN96  values, 
ratios,  and  final  average  ratio. 


line3 

Brault  & 

Brown68t> 

HI TRAN 9 6b 

ratio 

023P22 

2.44E-27 

2.65E-27 

0 . 92145 

019P18 

6. 47E-27 

6 . 87E-27 

0 . 94205 

P23P23 

5. 67E-27 

5 . 56E-27 

1.019234 

P23Q22 

6. 99E-27 

6.12E-27 

1 . 142344 

013P12 

1. 67E-26 

1.62E-26 

1.030864 

P21P21 

9. 10E-27 

9.36E-27 

0 . 972118 

P21Q20 

1 . 02E-26 

1 .04E-26 

0 . 982659 

P19P19 

1.43E-26 

1.47E-26 

0 . 970808 

O11P10 

1 . 90E-26 

1 . 78E-26 

1.065022 

P19Q18 

1 . 70E-26 

1 . 65E-26 

1.030303 

P17P17 

2.06E-26 

2 . 16E-26 

0.954588 

0  9P  8 

1 .75E-26 

1 .72E-26 

1.018034 

P15P15 

2 . 75E-26 

2 . 93E-26 

0 . 937287 

P15Q14 

3.45E-26 

3 . 38E-2  6 

1.019805 

average 

1.000469 

a)  Branch  symbol  for  AN,  N",  Branch  symbol  for  AJ,  J" 

b)  in  units  of  cnrl/(molec*cm'2) 
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Table  11. 


Measured  band  intensities,  derived  Einstein-A 
coefficients,  and  statistical  degeneracy  factors  for  the 
a’Ag  (v=0)  X3l-g  (v=0)  band. 


Reference 

Soo 

cm-Ulmolec'Ctrr2) 

Ae'AV->e"A"v" 

sec-1 

d//du 

Badger  et  al.65 

3.66  x  10-24 

2.58  x  lO'4 

3/2 

Lin  et  al.66 

9.4  x  10'24 

1.3  x  lO’4 

3/1 

Hsu  et  al.67 

2.1  x  lO'24 

2.9  x  lO*4 

3/1 

fit  of  data68  (t) 

3.69  x  lO'24 

2.59  x  lO*4 

3/2 

T  Adopted  for  HITRAN  96 


The  ai-Ajr  fv-n  <-  X3.£-£  (v-0)  band 

The  molecular  constants  for  the  a*Ag  (v=l)  state  are  from 
Brault.7^  Line  positions  have  changed  0.001302  cm-1  on  average. 
The  Einstein-A  coefficient  used  is  1/200  of  the  value  of  the  aJAg 
(v=0)  <-  X3£-g  (v=0)  band,76  hence  the  line  intensities  are  increased 
by  1.46  on  average.  Unfortunately,  to  our  knowledge,  no  other 
observations  of  this  band  are  available  to  perform  a  proper  test  of 
this  adopted  ratio.  The  line  positions  and  energies  have  only 
changed  by  -0.000273  cnr1  and  0.000388  cm*1  on  average. 

The  alA^  fv=01  XlX-£  (v=l)  band 

The  wavenumbers  have  changed  due  to  the  reformulation  of 
the  energy  expressions  for  both  the  upper  and  lower  states, 5  3, 6  3 
resulting  in  an  average  change  of  0.0157  cm*l.  The  line  intensities 
are  calculated  using  one-tenth  of  the  Einstein-A  coefficient  of  the 
a1  Ag  (v=0)  X3I-g  (v=0)  band76  and  are  larger  than  previous 


HITRAN  databases37  by  an  average  of  2.04.  The  same  comment  as 
above  applies  to  the  adopted  ratio  used  for  this  band. 


The  bl£±£L  (v  =01  <-  Xi£i£  (v=0)  band  (A-band) 

The  constants  for  the  Mz+g  (v=0)  state  are  from  Zare  et  al.77 
and  those  for  the  X3L*g  (v=0)  state  are  from  Rouille  et  al.53.  The 
wavenumbers  differ  by  -0.00127  cm*1  on  average,  with  the 
maximum  difference  being  0.00383  cm-'.  The  average  energy 

difference  is  0.00134  cm-1.  Intensity  ratios  vary  from  0.825  to 
0.877,  with  the  average  ratio  being  0.862.  The  Einstein-A  coefficient 
has  been  changed  from  the  value  of  Miller  et  al.,4^  0.077  sec-1  to  the 
newer  value  of  0.0887  from  Ritter  and  Wilkerson.71 

The  bl£%  (v=n  <-  X^I-g,.  (v=0)  band  (B-band) 

The  constants  for  the  b]I+g  (v=l)  state  are  from  Zare  et  al.77 
and  those  for  the  X3I'g  (v=0)  state  are  from  Rouille  et  al.53  The 
wavenumber  and  energies  arise  from  changes  made  to  the  ground 
state  energies.  The  maximum  wavenumber  and  energy  differences 
are  0.00383  cm'1  and  0.00370  cnr1,  respectively.  The  intensity 
ratios  range  from  0.782  to  0.822  with  an  average  of  0.811.  The 
value  of  the  Einstein-A  coefficient  is  from  Giver  et  al.,48  A(l- 
0)=0. 00724  sec-1,  but  was  incorrectly  used  in  the  previous  program 
where  it  was  set  to  A(l-0)=0. 00591  sec"1. 


The  bll±£  fv=2i  <-  Xl£-£,  fv=0)  band  (y-band) 

The  constants  for  the  b^+g  (v=2)  state  are  also  from  Zare  et 
al.77  and  those  for  the  X3E‘g  (v=0)  state  are  from  Rouille  et  al.53.  The 
average  wavenumber  and  energy  differences  are  0.00694  cm'1  and 
0.00104  cm-1.  The  intensity  ratios  range  from  1.22  to  1.27.  The 
Einstein-A  coefficient  of  Melieres  et  al.50  is  now  used,  A(2- 
0)=2.69xl0'4  sec1.  Compared  with  the  old  value  from  Miller  et  al.49 
of  A(2-0)=3.24xl0'4  sec1  (which  was  incorrectly  used  as  A(2- 
0)=2.2xl0'4  sec-1  in  the  HITRAN8237  calculations)  an  average  ratio  of 
0.813  is  obtained. 

The  blL±£  (v=ll  <-  X^L-g.  £vhU  band 

The  molecular  constants  used  in  the  calculation  of  energies  are 
from  Zare  et  al.77  for  the  b^+g  (v=l)  state  and  from  Rouille  et  al.53 
for  the  X3Z'g  (v=l)  state.  The  line  positions  show  a  maximum 
difference  of  0.0857  cm4  with  an  average  of  0.0207  cm'1.  The 
difference  is  due  to  the  change  in  the  ground  state  energies.  The 
Einstein-A  coefficient  is  from  Giver  et  al  48  The  intensity  ratios  range 

from  0.962  to  0.999.  These  changes  are  due  to  the  correct  partition 

sums  and  removal  of  some  approximations  in  the  formulas. 

The  blZ±E  (v=01  <-  XlZ-^  (v=l)  band 

The  molecular  constants  are  from  Zare  et  al.77  for  the  b^+g 
(v=0)  state  and  from  Rouille  et  al.53  for  the  X3Z'g  (v=l)  state.  The 

changes  in  the  line  positions  are  from  the  new  lower  state  energies. 

The  maximum  value  is  0.0333  cm*1  and  the  average  is  0.00935  cm*1. 
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The  Einstein-A  coefficient  is  from  Galkin.47  The  resulting  intensity 
ratios  range  from  0.968  to  1 .00. 


2.4.2  The  160180  species 

The  Xl£-&.  (v=01  <-  Xl£-t  (v=0)  band 

The  present  calculations  use  the  improved  molecular  constants 
of  Mizushima  and  Yamamoto.78  The  resulting  energies  differ  from 
previous  calculations  by  0.0174  cm4  on  average  with  a  maximum 
difference  of  0.136  cm-1  at  N=53.  The  average  difference  in  line 
position  is  0.00112  cm-1  with  the  largest  difference  being  0.00962 
cm-1.  The  line  intensities  differ  only  by  a  few  percent  maximum 
which  is  attributed  to  the  improved  partition  sums. 

The  alA£  (v =01  f-  fv=01  band 

There  is  a  slight  change  in  some  of  the  wavenumbers  due  to 
the  change  in  energy  formulation  for  the  X^L’g  (v=0)  state. 53  This 
gives  rise  to  a  maximum  difference  of  0.00580  cm"1  and  an  average 
difference  of  -0.00127  cm-1.  Intensity  ratios  are  ~0.5  as  expected 
and  49  of  the  previous  lines  do  not  make  the  intensity  cut  off. 

The  b-L£— (t  (v=0)  j —  fv— 0)  band 

The  energies  of  the  X3Z"g  (v=0)  state  are  calculated  using  the 
constants  of  Mizushima  and  Yamamoto.78  The  constants  for  the 
b^+g  (v=0)  state  are  from  Babcock  and  Herzberg.79  Wavenumber 
differences  of  0.0242  cnr1  at  N"=34  are  observed  with  an  average 


difference  of  -0.00577  cm-1.  The  average  intensity  ratio  is  0.882  due 
mostly  to  the  change  in  partition  functions. 


The  blS+^r  (v=n  <-  Xll-£,  (v=01  band 

The  constants  for  the  t^L+g  (v=l)  state  are  from  Benedict.80 
Those  for  the  X3Z-g  (v=0)  state  are  from  Mizushima  and  Yamamoto.78 
The  maximum  wavenumber  differences  is  0.0107 
cm-1.  This  difference  is  due  to  the  change  in  the  lower  state 
energies.  The  Einstein-A  coefficient  is  from  Giver  et  al.48  The 
average  ratio  of  the  intensities  is  0.830. 

The  b— S— (v=21  j —  ( v— 0)  band 

The  constants  for  the  b^+g  (v=2)  state  are  from  Zare  et  al.77 
Those  for  the  X3Z-g  (v=0)  state  are  from  Mizushima  and  Yamamoto.78 
The  maximum  wavenumber  difference  is  0.150  cm*1  and  the  average 
intensity  ratio  is  0.837.  Caution  must  be  used  in  interpreting  these 
numbers  since  the  comparison  is  based  on  only  3  lines. 

2.4.3  The  160170  species 


The  XlS-a  (v=0)  <-  XlS-£  (v=0)  band 

The  data  for  this  band  are  taken  from  the  Jet  Propulsion 
Laboratory  catalogue.81  There  are  10,787  lines  in  the  JPL  file.  The 
data  were  filtered  through  the  wavenumber  dependent  cutoff 
resulting  in  2601  lines  from  0.000012  cm-1  to  186.15  cnr1  in  the 
final  file.  Note  the  isotopic  abundance  factor  was  inadvertently 
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omitted  from  the  1996  HITRAN  database.  Thus,  the  ratio  of  the  line 
intensity  (S92/S96)  is  0.000750  on  average.  In  order  to  properly 
use  the  intensities  for  these  data,  they  should  be  multiplied  by  Ia  = 
0.000742235.  The  line  positions  and  lower  state  energies  have  only 
changed  slightly,  0.000041  cm-1  and  -0.000352  cm-1  average 
difference  respectively. 

The  bl£±g  (v=ll  <-  Xll-jr  (v=0)  band 

These  data  are  from  Benedict  and  Brault82  and  have  not 
changed  from  the  1982  HITRAN  database. 


2.5  Other  Changes 

We  have  added  reference  and  error  codes  to  the  line  parameter 
database.  The  error  code  (see  HITRAN96  manual83)  for  the 
halfwidths  is  set  to  4.  The  error  code  for  the  line  positions  of  X3L'g 
<-  X3Z'g  electronic  band  is  set  to  4  and  all  other  error  codes  are  not 
utilized  on  the  database  (i.e.  set  to  0).  We  have  also  labeled  the 
electric  quadrupole  and  magnetic  dipole  transitions  by  the  lower  case 
letters  q  and  d,  respectively,  in  the  sym  field  of  the  rotational 
quantum  number  character  string,  i.e.  Br,  F",  Br,  N",  Br,  J", _ ,  Sym. 

These  data  are  available  in  the  1996  HITRAN  database.83 

2.6  O2  Continuum  Absorption 

It  is  known  that  the  02  X3I-g(v")  -  a^gfv')  absorption  bands 
exhibit  both  discrete  (rotational)  line  structure  and  pressure-induced 


continuous  absorption. 65  The  most  important  bands  are  the  v"-0, 
v’=0  at  1.27  pm  (7882  cm'1),  v"=0,  v'=l  at  1.06  pm  (9366  cm*1),  and 
v"=l,  v'=0  at  1.6  pm  (6326  cm'1).  While  the  rotational  lines  of  the 
(0-1)  band  are  weaker  than  those  of  (0-0)  band,  and  those  of  (1-0) 
band  are  weaker  than  the  (0-1)  band,  the  continuum  absorptions  are 
of  more  similar  intensity. 

During  the  update  of  the  (0-0)  band  line  parameters, 
theoretical  calculations  were  compared  with  absolute  atmospheric 
transmittance  obtained  with  the  University  of  Denver  Absolute  Solar 
Transmittance  Interferometer  (ASTI).84  The  results  show  good 
agreement  of  the  line  structure  but  a  clear  indication  of  the 
underlying  continuum  (not  modeled),  with  the  P,  R  (no  Q)  shape  of 
the  envelope  under  the  absorption  lines.  The  continuum  is  clearer  at 
the  higher  spectral  resolution  (>0.5  cm  OPD). 

More  recent  ASTI  data  in  the  9400  cm-1  and  6400  cm'  1 
regions  also  show  a  strong  continuum,  similar  to  that  described 
above.  It  is  thus  proposed  that  this  is  due  to  the  pressure-induced 
absorption  of  the  v"=0,  v’=l  and  v"=l,  v’=0  bands  respectively.8 5 
Laboratory  data86  are  consistent  with  these  conclusions,  but  indicate 
no  pressure-induced  absorptions  under  the  X3Z'g(v")  -blL+g(v  ) 
bands.  The  pressure-induced  absorption  in  the 

X3£-g(v"=0)  -  X3Z~g(v'=l )  has  been  well  documented  in  previous 

publications.59*87’88 

In  the  HITRAN  database  we  have  included  the  individual  line 
parameters  of  these  bands,  but  no  cross-sections  are  provided  for  the 
continuum.  These  cross-sections  will  be  forthcoming  on  the  next 
version  of  the  HITRAN  database. 


3.0  HALFWIDTHS  AND  LINE  SHIFTS  OF  H20  TRANSITIONS 


Water  vapor,  the  most  important  greenhouse  gas  in  the 
terrestrial  atmosphere,  is  the  principle  atmospheric  absorber  of 
infrared  radiation  with  some  50,000  spectral  transitions  ranging 
from  the  microwave  to  the  visible  portions  of  the  spectrum. 3  In 
order  to  interpret  remote  sensing  measurements  and  understand 
radiative  properties  of  the  atmosphere,  it  is  useful  to  know  spectral 
parameters  along  with  their  temperature  dependence.  For  the 
terrestrial  atmosphere,  the  collisional  interaction  of  water  vapor  with 
nitrogen  and  oxygen  must  be  understood. 

A  knowledge  of  the  fundamental  spectroscopic  parameters  for 
H20  in  the  presence  of  other  gases  (generally  referred  to  as  the 
perturbing  gas  or  buffer  gas)  is  important  for  applications  where 
these  gases  are  the  principal  components  of  the  atmosphere  or 
system  in  question.  For  example,  the  Earth’s  nearest  neighbors  - 
Venus  and  Mars  -  have  atmospheres  composed  predominantly  (at 
the  95%  level)  of  carbon  dioxide.  Spectroscopic  observations  of  these 
atmospheres  from  spacecraft  and  ground-based  platforms  can 
provide  much  information  about  the  presence  and  abundance  of 
water  vapor.  Water  vapor  and  carbon  dioxide  are  also  the  principal 
products  of  combustion,  thus  spectroscopic  techniques  can  be  used  to 
monitor  combustion  processes. 

The  atmospheres  of  the  giant  planets  are  made  mostly  of  H2 
and  He,  the  same  constituents  that  are  the  dominant  constituents  of 


the  Sun.  This  composition  reflects  the  ability  of  the  massive  giant 
planets  to  have  gravitationally  captured  large  amounts  of  gases  from 
their  birthplace  in  the  primordial  solar  nebula  as  well  as  solids  (e.g., 
water)  that  were  eventually  volatilized  into  heavier  gases  in  their 
atmospheres.  By  obtaining  accurate  estimates  of  the  water 
abundance  in  the  atmospheres  of  the  giant  planets,  constraints  can 

be  placed  on  the  relative  efficiency  with  which  they  accumulated 

gases  and  solids  from  the  primordial  solar  nebula.  Thus,  for  these 

applications  it  is  important  to  know  the  spectral  parameters  for 

water  interacting  with  hydrogen  and  helium. 

The  fundamental  spectral  parameters  needed  for  a  retrieval 

are  the  line  position,  intensity,  collision-broadened  halfwidth, 
pressure-induced  line  shift,  and  lower-state  energy. 3  The  least  well 

characterized  of  these  parameters  are  the  pressure  induced 

halfwidths  and  line  shifts  for  the  collision  partners  mentioned  above. 
The  need  of  the  spectroscopic  and  remote  sensing  communities  is  to 
know  such  parameters  to  an  accuracy  of  l-5%.89'90 

Theory  can  contribute  to  this  purpose  in  two  ways.  First,  the 
vast  spectral  range  and  environmental  conditions  are  difficult  to 
cover  exhaustively  in  the  laboratory.  Second,  nonlinear  pressure- 

broadening  effects  require  a  solid  theoretical  model  in  order  to 
unravel  the  sometimes  complicated  spectra  observed.  Here  we  focus 
on  improving  the  basic  theoretical  apparatus  for  calculating 
halfwidths  and  line  shifts  for  larger  species  present  in  the  terrestrial 
atmosphere,  e.g.  H2O,  O3,  CH4,  etc.  The  model  chosen  is  a  complex 
implementation  of  the  semiclassical  formalism  of  Robert  and 
Bonamy91  (RB)  for  several  reasons. 
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1)  The  formalism  is  based  on  the  resolvent  operator  formalism  of 
Baranger,92  Kolb,93  and  Greim94  (BKG).  The  application  of  linked- 
cluster  techniques95  to  the  BKG  formalism  leads  to 
developments91  >96-98  which  eliminate  the  awkward  cutoff  procedure 
that  characterized  earlier  theories. 99-102  Also,  the  cumulant9  3 
expansion,  as  incorporated  in  the  RB  formalism,  provides  a  better 
description  of  the  long  term  dynamics  of  the  collisions. 

2)  The  formalism  is  complex  valued,  yielding  halfwidths  and  line 

shifts  from  a  single  calculation. 

3)  The  intermolecular  dynamics  are  treated  more  realistically  than 

in  earlier  theories,  i.e.  using  curved  rather  than  straight  line 
trajectories.  This  has  important  consequences  in  the  description  of 
close  intermolecular  collisions  (small  impact  parameters). 

4)  Connected  to  item  (3)  is  the  incorporation  in  the  RB  theory  of  a 

short  range  (Lennard-Jones  6- 1 21 03 )  atom-atom  component  to  the 
intermolecular  potential.  This  component  has  been  shown  to  be 
essential  for  a  proper  description  of  pressure  broadening,  especially 
in  systems  where  electrostatic  interactions  are  weak.104  (Here,  the 
notion  of  strong  and  weak  collisions  adopts  the  definition  of  Oka.105) 

5)  This  formalism  allows  the  removal  of  all  "adjustable"  parameters 

so  as  to  arrive  at  a  more  predictive  theory. 

While  not  as  rigorous,  this  formalism  has  advantages  over  more 
exact  theories106-111  in  that  it  is  tractable  for  molecules  of 
atmospheric  interest  such  as  H2O,  O3  and  CH4.  The  emphasis  has 
been  to  develop  a  formalism  which  can  consider  such  atmospheric 
species  and  produce  results  at  the  5%  uncertainty  level  as  required 
by  the  spectroscopic  and  remote  sensing  communities. 


The  new  feature  of  this  work  is  the  development  of  the 

complex  Robert-Bonamy  formalism  (CRBF)  which  can  be  applied  to 
larger  atmospheric  species.  Previously,  the  imaginary  components  of 
the  resonance  functions  that  arise  in  the  theory  were  only 

determined  for  a  small  number  of  terms  that  occur  in  the 
electrostatic  potential.112*113  We  have  determined  these  resonance 
functions  for  all  integrals  that  arise  in  a  general  spherical  tensor 
expansion  of  the  potential.114*11^ 

It  should  be  noted  that  within  the  complex  Robert-Bonamy 

formalism,  the  imaginary  parts  of  the  S  matrix  expansion  affect  both 
the  calculation  of  the  line  shift  and  the  halfwidth,  an  effect  not 

achieved  in  Anderson-Tsao-Curnutte  (ATC)  theory. 99-Hn  Calculations 
based  on  the  complex  Robert-Bonamy  formalism  indicate  that  the 
effect  of  the  imaginary  terms  is  very  important  in  determining  the 
halfwidth. 1 14-1 1 8  The  effect  of  the  imaginary  components  on  the 
halfwidths  varies  from  transition-to-transition  and  perturber-to- 
perturber  but  can  be  as  much  as  25%.  The  change  is  generally 
(almost  always)  in  the  direction  of  better  agreement  with 
experiment. 

In  this  work,  we  summarize  our  CRBF  calculations  of  H2O 
pressure-broadened  halfwidths  and  pressure  induced  line  shifts  by 
several  perturbing  gases:  N2,  O2,  CO2,  and  H2,  and  compare  the  results 
with  measurements.  We  have  attempted  to  compare  with 

measurements  where  both  halfwidths  and  line  shifts  were  made 
simultaneously,  although  this  was  not  always  possible  for  the 
perturbing  gases  listed  above.  For  N2  and  O2  as  the  buffer  gases,  the 
data  are  those  of  Grossmann  and  Browell119  for  the  3vi  +  V3  and 


2vi+2v2+V3  bands  of  H2O.  There  are  several  reasons  for  this  choice: 
the  halfwidths  from  this  work  exhibited  high  quality  in  a  survey  of 
measurements  of  H2O  halfwidths,120  and  the  line  shifts  for  these  ro- 
vibrational  transitions  are  sufficiently  large  in  magnitude  to  be 
measured  with  reasonable  certainty.  Measurements  of  the  halfwidth 
and  line  shift  of  transitions  in  the  Vj,  2V2,  and  V3  bands  of  H2O  in  a 
CO2  buffer  gas  are  from  Gamache  et  al.116-121  Measurements  of  630 
room-temperature  H2-broadened  halfwidths  of  H2O  transitions  in  the 
pure  rotational,  and  three  fundamental  bands  between  55  and  4045 
cm*1  are  those  reported  by  Brown  and  Plymate.122  Together,  these 
data  amount  to  772  transitions  for  which  halfwidth  data  are 
available  and  212  transitions  for  which  line  shift  data  are  available. 
The  CRBF  calculations  are  made  at  a  temperature  of  296  K  and  are 
discussed  below.  Note,  because  all  measurements  are  within  a  few 
degrees  of  296K  no  temperature  correction  is  applied  to  the  data. 

3.1  Complex  Robert-Bonamy  Formalism 

Within  the  Robert-Bonamy  formalism,  the  line  shift,  5,  and 
halfwidth,  y,  for  the  transition  f<-  i  are  given  by  the  real  and  minus 
the  imaginary  parts,  respectively,  of  the  diagonal  elements  of  the 
complex  relaxation  matrix.123  In  computational  form,  the  line  shift 
and  halfwidth  are  usually  expressed  in  terms  of  the  Liouville 
scattering  matrix,92-124 

6f<_i  £  <J2lp2U2>f  2rcb  db  sin(Si+/  m{S2})  cRe^ 

2%c  j2  J0 


(22) 


Yf«_j  =  —  V-  X  <J2lp2IJ2)f  2nb[l  -  cos(Si+/ w(S2j)  (23) 

2*C  h  Jo 

where  v  is  the  mean  relative  thermal  velocity,  p2  and  n2  are  the 
density  operator  and  number  density  of  perturbers,  and  b  is  the 
impact  parameter.  Si  and  S2  are  the  first-  and  second-order  terms  in 
successive  approximation  of  the  Liouville  scattering  matrix  and 

depend  on  the  intermolecular  potential. 

The  corresponding  equations  from  Anderson-Tsao-Curnutte 

theory99"1^1  are 

5f*_i  =  -^2JL  £  <J2lp2IJ2>f  2rcb[Si+/m(S2}]  db  (24) 

2rcc  j2  A 

Yf^.i  X  <J2lp2IJ2>|  2nb  i?e(S2}  db  (25) 

2nc  h  Jo 


where  the  variables  are  as  in  Eqs.  (22)  and  (23).  Note,  in  ATC  theory, 
the  halfwidth  is  the  integral  over  the  components  of  the  S2  term,  or 
equivalently,  the  sum  of  the  integral  of  each  term.  (In  general,  for 
each  term  in  the  potential  the  S2  term  is  composed  of  3  operators, 
see  Ref.  100  and  101  for  details.)  The  similar  feature  of  additivity 
also  holds  for  the  line  shift  in  ATC  theory,  Eq.(24).  This  is  not  true  in 
the  theory  of  Robert  and  Bonamy  because  the  method  of  cumulants95 
places  the  Liouville  scattering  matrix  in  the  exponential.  More 
important,  when  the  real  and  imaginary  components  of  the 
relaxation  operator  are  taken,  one  arrives  at  Eqs.  (22)  and  (23)  in 
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which  the  real  and  imaginary  components  of  the  expanded  S  matrix 
are  present  in  both  expressions.  All  previous  RB  calculations121 -125' 
132  for  asymmetric  rotor  molecules  employing  the  atom-atom  and 
electrostatic  potentials  have  assumed  the  imaginary  contribution  to 
the  halfwidth  to  be  negligible.  The  effects  of  the  imaginary 
components,  the  cosine  term  in  Eq.  (23),  on  the  halfwidths  for  the 
perturbing  gases  considered  here  are  summarized  below. 

In  the  CRBF  formalism  the  intermolecular  potential  is  a 
combination  of  three  kinds  of  terms:  electrostatic,  atom-atom,  and 
induction-dispersion.  The  first  of  these  is  based  on  the  permanent 
distribution  of  charge  characterizing  a  given  molecule,  usually 
expressed  as  an  expansion  in  moments  (dipole,  quadrupole,  etc.). 
Gray  and  Gubbins133*134  have  shown  that  the  interaction  resulting 
from  two  such  charge  distributions  in  proximity  can  be  expressed  in 
the  form  of  a  spherical  tensor  expansion,  presented  here  in  a  slightly 
modified  form  for  a  H20-linear  system, 

v  =  y  y  Y  niw(i) 

nl  w,q  j^q+^]+^2+2w 
i  mj  m2 
m 

(26) 

®  CUi  £2  £\  mi  m2  m)D^,ini(Qi)  D^20(Q2)  Yto(co) 

where  C(h  12  l;mi  m2  m)  is  a  Clebsch-Gordan  coefficient,  Qj=(ai,  Pi,  Yi) 
and  Q2=(a2,  p2,  y2)  are  the  Euler  angles  describing  the  molecular  fixed 
axis  relative  to  the  space  fixed  axis.  Subscripts  1  and  2  refer  to  the 
radiating  and  perturbing  molecules,  respectively,  and  to  =  (0,(j)) 


describes  the  relative  orientation  of  the  centers  of  mass  (see  Figure 
9).  As  stated  previously,  R  is  the  center  of  mass  separation.  The 
coefficients  U(...)  depend  on  electrostatic  moments  of  the  radiator 
(H2O)  and  perturber  gases  according  to  expressions  given  in  appendix 
A  of  ref.  115. 

Although  an  electrostatic  potential  is  suitable  for  large 
intermolecular  separations,  it  is  clearly  inadequate  at  short  range  for 
neutral  molecules.  For  example,  the  electrostatic  terms  in  Eq.  (26) 
are  always  orientation-dependent,  whereas  molecules  in  any 
orientation  must  ultimately  repel  one  another  at  close  range.  The 
required  terms  are  isotropic,  i.e.  with  1 1  =  1 2  =0,  which  cannot  be 
supplied  by  electrostatics. 


Figure  9  Geometry  for  H2O-A2  system.  R  is  the  center  of  mass 
separation  and  o)  is  the  relative  orientation.  The 
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molecular  fixed  axes  are  defined  relative  to  the  space 
fixed  axis  by  Euler  angles. 

To  correct  this  deficiency,  an  atom-atom  potential  is  added  to 
the  total  interaction.  This  potential  is  defined  as  the  sum  of  pair¬ 
wise  Lennard-Jones  6-12  interactions  between  atoms  of  molecules  1 
and  2, 
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The  subscripts  li  and  2j  refer  to  the  i th  atom  of  molecule  1  and  the 
)th  atom  of  molecule  2,  n  and  m  are  the  number  of  atoms  in 
molecules  1  and  2,  respectively,  and  ey  and  Oij  are  the  Lennard-Jones 
parameters  for  the  atomic  pairs.  The  heteronuclear  atom-atom 
parameters  ey  and  Gy  are  then  constructed  from  homonuclear  atom- 
atom  parameters  (et  and  Gj)  by  the  "combination  rules"135 
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Gray  and  Gubbins133  have  shown  how  the  atom-atom  potential 
may  be  expressed  in  the  form  of  Eq.  (26),  and  provide  explicit 
expressions  for  U(...)  in  terms  of  £y  and  Gy.  With  such  expressions, 
CRBF  theory  can  be  made  to  incorporate  any  combination  of 
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electrostatic  and  atom-atom  components  one  wishes  to  employ.  The 
natural  choice  is  simply  to  add  the  two,  i.e., 

V(total)  =  V(elec)  +  V(atom-atom)  (3.8) 

using  available  literature  values  for  all  electrostatic  and  atom-atom 
parameters.  It  has  been  shownl 1 7*1 1 8”1 36  that  the  addition  of  the 
atom-atom  potential  to  the  overall  potential  does  not  represent  any 
"over  counting"  or  duplication  of  potential  interactions. 

Before  taking  up  discussion  of  the  third  interaction  (induction- 
dispersion),  three  further  points  should  be  emphasized  with  regard 
to  the  atom-atom  potential. 

1.  When  expressed  in  the  form  of  Eq.  (26),  the  atom-atom 
potential  can  be  understood  as  two  simultaneous  expansions.  One  is 
defined  by  the  tensorial  ranks  1]  and  h,  which  determine  the 
symmetry  of  the  interaction.  A  second  expansion  is  defined  by  the 
sum  li+l2+2w,  which  we  call  the  order  of  the  expansion  because  it  is 
the  integer  added  to  q  in  Eq.  (26).  It  hardly  needs  to  be  emphasized 
that  any  calculation  of  halfwidth  or  line  shift  should  be  converged 
with  respect  to  both  order  and  rank. 

2.  The  atom-atom  potential  contains  an  isotropic  part  (1 1  =  1 2=0)  which 
in  general  are  comprised  of  terms  of  the  form  1/R6,  1/R8,  1/R10,  etc. 
In  our  implementation,  these  are  used  to  define  the  intermolecular 
trajectory.  For  the  purpose  of  comparison,  however,  it  is  convenient 
to  fit  this  summation  to  an  effective  isotropic  (or  "hetero-molecular") 
Lennard-Jones  6-12  potential. 
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(30) 


y  iS0 


The  approximation  represented  by  Eq.  (30)  had  no  significant 
impact  on  our  calculations  other  than  to  simplify  the  trajectory 
calculations  somewhat. 

3.  The  remaining  parts  of  the  atom-atom  potential  are  anisotropic 
(non  zero  li  or  I2).  These  terms  may  be  discussed  with  the  same 
lexicon  that  one  applies  to  electrostatic  interactions.  For  example, 
when  1 1  =1  and  12=2,  we  refer  to  a  "dipole-quadrupole"  interaction.  It 
should  be  kept  in  mind  that  these  are  symmetry  appellations  only, 
and  do  not  refer  to  electrostatics.  Thus,  it  happens  that  between  H2O 
and  N2,  for  example,  the  lowest  anisotropic  atom-atom  term  is 
"dipole-monopole",  i.e.  with  1]  =  1  and  12=0.  Such  terms  can  contribute 
significantly  to  the  line  broadening  mechanism.1 17*' 18*136 

In  a  third  category  of  interactions  relevant  to  intermolecular 
collisions  are  the  isotropic  induction  and  dispersion  potentials, 
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(31) 


yd  i  s  _  2_  F] E2  ^1^2 

iso  *  2  E,+E2  R6 


(32) 


where  p.]  is  the  dipole  moment  of  water  vapor  and  a  and  E  are  the 
polarizability  and  ionization  potentials  for  water  vapor  and  its 
perturber.  By  design,  these  interactions  occur  with  the  same  power 
of  1/R,  and  in  some  sense  represent  the  same  forces,  as  the  attractive 
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part  of  Eq.  (30).  Thus,  they  play  no  role  in  our  trajectory 
calculations,  although  they  do  serve  as  a  convenient  check:  we  find 
that  the  magnitude  of  the  combined  induction-dispersion  interactions 
is  typically  about  one-half  the  value  of  the  corresponding  atom-atom 
term.  This  is  consistent  with  observations  of  Gray  and  Gubbins,133 
who  point  out  that  Eqs.  (31)  and  (32)  are  themselves  merely  the  first 
terms  in  an  expansion,  and  by  Lynch  et  al.115.  In  principle,  one  could 
arrive  at  the  vibrational  dependence  of  the  intermolecular  atom- 
atom  potential  using  contact  transformations  of  Eq.  (27)  to  the 
appropriate  vibrational  state.  However,  the  isotropic  induction  and 
dispersion  potentials  provide  a  convenient,  experimentally 
constrained  method  of  computing  the  difference  in  isotropic  potential 
for  different  vibrational  states,  via  the  vibrational  dependence  of  p.] 
and  a\  (this  difference  appears  in  the  Si  term  of  CRBF  theory,  as 
described  below).  This  difference  is  not  accounted  for  in  any  other 
part  of  the  theory,  and  represents  no  duplication  of  potential 
interactions. 

The  first  order  term  in  Eqs.  (22)  and  (23)  depends  only  on  the 
difference  in  the  isotropic  part  of  the  interaction  potential  between 
the  initial  and  final  vibrational  states,  v;  and  vf,  of  the  radiator 


Si 


iso 


lVf  >  -  <Vjl  Viso  I v j  >] 


(33) 


where  the  time  integral  is  evaluated  along  the  relative  trajectories  of 
the  radiator  and  perturber.  The  vibrational  dependence  of  the 
isotropic  interaction  is  approximated  from  the  isotropic  induction  and 
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dispersion  potentials,  Eqs.  (31)  and  (32).  The  dependence  of  the 
dipole  moment  and  polarizability  on  the  vibrational  state  of  the 
molecule137’ *38  is  considered  in  the  calculations.  For  a  ro-vibrational 
transition  of  the  radiating  molecule,  the  first  order  term  can  be 
written, 


Si  = 

Rc  vc  8 

where  Rc  and  vc'  are  the  distance  and  effective  velocity  at  the  point 
of  closest  approach  within  the  parabolic  approximation.91 

The  vibrational  dependence  of  the  dipole  moment  of  water 
vapor  was  investigated  by  Shostak  and  Muenter139  and  is  given  in 
Debyes  by 


~-~~-(al,f  ’  al,i  )  '  M.f  '  ^l.i) 
2  / 1 +/2 


(34) 


p  =  1.855  +  0.0051v,  -  0.0317v2  +  0.0225v3  ,  (35) 

where  v;  is  the  number  of  quanta  in  the  ith  normal  mode.  The 
vibrational  dependence  of  the  polarizability  of  water  vapor  was 
obtained  by  Luo  et  al.140  and  is,  in  atomic  units,141 

a  =  9.86  +  0.29v,  +  0.03v2  +  0.28v3  .  (36) 

The  second-order  term  is  the  complex  analog  of  that  appearing 
in  the  familiar  ATC  theory,99'101 

*  c\ 

S2  =  S2J2  +  $2/2  +  S2, middle 


6  1 


(37) 


where  the  notation  is  that  of  Anderson.100  The  case  in  which  the 
imaginary  part  of  S2  is  ignored  has  been  discussed  thoroughly  in  the 
literature. "-,01>142>143  Note,  S2, middle  only  has  a  real  component. 
The  other  terms  are  complex  functions  and  can  be  written  in  the 
form 

S2,f2= - - -  X  X  D(^2,nanb,JfJfJ2J2)  F^l2nb  (cof2>f2')  (38) 

H2  [Jf]  [J2]  V?  nanb 
J2Jf 

where  [J]=2J+1,  n  =  (ni,n2)  and  C0f2,r2’  =  (Ep  -  Ef  +  E2’  -  E2)  where  Es  is 
the  energy  of  the  state  f,  f,  2’,  or  2.  S2,i2  is  obtained  from  Eq.  (3.17) 

by  replacing  f  with  i.  The  D  terms  are  reduced  matrix  elements  for 
the  internal  states  of  the  radiator  and  perturber.1 15  The  F  terms  are 
the  resonance  functions  and  are  given  by 

F  nf2nb  (co)  =  -  *  1  X  X  C(^i^a,m1m2ma)C(^,^b,m1m2mb) 

KlJ  l*'2j  w»  wb 

qa  qb  ma  mb 
m  1  m2 

UUi^a,  na,waqa)  U*(^i^2^b,  nb,wbqb)  (39) 

im,  Yw[m(t )]  T  dt,  e.imf  Y*bmb[co(t')] 

Rp*(t)  L  Rpb(t’) 

where  pa=qa+^i+^2+2wa  (similarly  for  pb)  and  the  other  terms  are  as 
described  earlier.  The  integration  is  along  the  relative  semiclassical 
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trajectory  of  the  molecules  involved.  All  integrals  are  evaluated,  real 
and  imaginary  terms, for  the  potential  discussed  above. 

3.2  Calculations 

Calculations  of  the  halfwidth  and  line  shift  based  on  the  CRBF 
method  are  made  for  transitions  of  H2O  perturbed  by  the  gases 

mentioned  above.  To  investigate  the  contribution  of  the  imaginary 

components  on  the  halfwidth  and  line  shift,  calculations  are  made 

with  the  different  imaginary  terms  and  compared.  The  different 
calculations  are  labeled  as  follows;  RRB:  the  real  components  of  the 
RB  formalism,  RS2:  the  real  components  and  the  imaginary  S2  terms, 
RSi:  the  real  components  and  the  Si  terms,  and  finally,  CRBF:  all  real 
and  imaginary  terms.  In  addition,  the  order  of  the  expansion 
(ll+l2+w)  of  the  atom-atom  potential  is  varied.  All  calculations  are 
made  at  the  HITRAN3  reference  temperature,  296  K.  No  temperature 
correction  is  applied  to  the  data  since  the  small  difference  in  the 
measurement  and  calculation  temperatures  is  negligible.  The  CRBF 
formalism  yields  the  halfwidth  and  line  shift  from  a  single 

calculation  and  there  are  no  adjustable  parameters  or  cutoff 
procedure  in  the  formalism. 

All  molecular  parameters  used  in  this  work  are  the  best 
available  values  from  the  literature.  No  molecular  constants  are 
adjusted  to  give  better  agreement  with  experiment.  The  dipole  and 
quadrupole  moments  of  water  vapor  are  taken  from  Refs.  139  and 
144,  respectively.  The  quadrupole  moment  of  nitrogen  is  from 
Mulder  et  al.145  and  that  for  oxygen  is  from  Stogryn  and  Stogryn146 
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There  have  been  a  number  of  measurements  of  the  quadrupole 
moment  of  carbon  dioxide147  which  range  from  4.0  to  roughly  4.6  in 
units  of  10'26  esu.  The  most  recent  measurement148  reports  a  value 
of  4.02(10)  x  10'26  esu  which  is  adopted  in  this  work.  For  hydrogen 
the  value  of  the  quadrupole  moment  is  from  Poll  and  Wolniewicz.149 
The  numerical  values  are  listed  in  Table  12.  The  vibrational 
dependence  of  the  dipole  moment  of  H2O  (Eq.(35))  is  taken  from 
Shostak  and  Muenter139  and  the  vibrational  dependence  of  the 
polarizability  of  H2O  (Eq.(36))  is  taken  from  the  work  of  Luo  et  al.14®. 
The  atom-atom  parameters.  Table  13,  were  obtained  using  the 
combination  rules,1 3^  (see  Eq.  28)  and  the  atom-atom  parameters  for 
homonuclear  diatomics  determined  by  Bouanich150  by  fitting  to 
second  virial  coefficient  data. 

Table  12.  Values  of  electrostatic  moments  for  the  water  vapor, 

N2,  O2,  CO2,  H2. 


Molecule 

Multipole  Moment 

Reference 

H20 

p  =  1.8549  x  10-18  esu 

139 

Qxx  =  -0.13  x  10'26  esu 

139 

Qyy  =  -2.5  x  10'26  esu 

139 

Qzz  =  2.63  x  10‘26  esu 

139 

n2 

Qzz  =  -1.4  x  10‘26  esu 

145 

O2 

0/./  =  -0.4  x  lO'26  esu 

146 

CO2 

Qzz  =  -4.02  x  10-26  csu 

148 

H2 

Qzz  =  -0.616  x  10-26  esu 

149 
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Table  13.  Values  of  the  heteronuclear  atom-atom  Lennard-Jones 
(6-12)  parameters  for  the  collision  pairs  considered  in 
this  work. 


Atomic  pair 

a/Angstrom 

e/kB  /  °K 

N-0 

3.148 

43.90  ! 

N-H 

2.990 

20.46 

0-0 

3.010 

51.73 

0-H 

2.850 

24.13 

H-C 

2.810 

32.18 

H-H 

2.683 

1 1.25 

0-C 

3.285 

40.43 

In  the  parabolic  approximation,  the  isotropic  part  of  the 
interaction  potential  is  taken  into  account  in  determining  the 
distance,  effective  velocity,  and  force  at  closest  approach.91  To 
simplify  the  trajectory  calculations,  the  isotropic  part  of  the  atom- 
atom  expansion  is  fit  to  an  isotropic  Lennard-Jones  6-12  potential 
(see  Eq.  30)  and  the  resulting  parameters  are  given  in  Table  14  for 
the  particular  collision  pairs  considered  in  this  work. 


Table  14.  Values  of  the  isotropic  atom-atom  parameters  used  to 
define  the  trajectories  of  the  collision  systems 
considered  in  this  work. 


Collision  System 

Order 

a/Angstrom 

e/kB  (°K) 

3.62 

96.6 

KWH''  -  » 

3.53 

105.22 

1 

Klin 

3.60 

164.35 

3.838 

119.3 

To  evaluate  the  reduced  matrix  elements  in  Eq.  (38),  ro- 
vibrational  wave  functions  are  needed.  For  the  3vi+V3  and  the 
2v  i  +  2v 2  + V 3  bands  the  constants  of  Grossmann  et  al.151  are  used. 
Because  of  limitations  in  computing  wavefunctions  for  high  J  states 
for  the  triad,  Vi,  2V2,  and  V3,  ground  vibrational  state  wavefunctions 
are  used  in  the  calculations  with  the  Watson-Hamiltonian  constants 
of  Flaud  and  Camy-Peyret.152  This  approximation  introduces  less 
than  a  percent  error  to  the  halfwidth  for  these  bands.  For  the 
rotation  and  V2  bands,  the  constants  of  Flaud  and  Camy-Peyret152  are 
used.  The  molecular  constants  for  N2,  O2,  and  H2  are  from  Huber  and 
Herzberg.153  For  carbon  dioxide,  the  molecular  constants  of  Rothman 
et  al.154  are  used. 

Finally,  the  ionization  potential  of  water  was  taken  to  be  a 
vibration  ally -in  dependent  12.6  eV.155  The  polarizability  and 

ionization  potential  of  the  perturbing  gases  are  given  in  Table  15. 

Table  15.  Values  of  the  polarizability  and  ionization  potential  for 
N2,  O2,  CO2,  H2. 


Perturbing 

gas 

Polarizability 
/10'24  cm3 

Reference 

Ionization 
potential  /ev 

Reference 

n2 

1 .740 

15.576 

157 

O2 

1 .580 

1  5  6 

12.063±0.001 

1  5  5 

CO2 

2.913 

156 

13.769±0.03 

158 

h2 

0.806 

156 

15.427±0.002 

159 
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3.3  Results  and  Discussion 


In  this  section,  we  present  a  summary  of  the  comparison  of  the 
calculated  and  measured  halfwidths  and  line  shifts  as  well  as 
describe  some  of  the  general  trends  observed.  There  are  several 
important  features  that  stand  out  from  this  work.  First,  we  address 
the  expansion  of  the  atom-atom  potential  and  the  effects  the  order  of 
expansion  has  on  both  the  halfwidth  and  line  shift.  The  order  of  the 
expansion,  defined  above,  is  varied  from  0  to  8  and  the  results 
compared.  Note,  the  electrostatic  and  induction-dispersion  potentials 
are  used  in  all  calculations  and  are  defined  as  the  0th  order 
calculation,  i.e.  no  atom-atom  terms.  Figure  10  shows  the  results  of 
calculations  of  the  halfwidth  for  transitions  belonging  to  the  3Vi+V3 
band  of  H2O  for  oxygen  as  a  perturbing  species.  Plotted  are  the 
halfwidth  versus  an  energy  ordered  index,  J"(J"+l)+Ka"-Kc"+l.  These 
include  the  measurements  of  Grossmann  and  Browell11^  with  error 
bars,  the  calculations  for  a  0lh  order  (diamond  symbols),  4th  order 
(circle  symbols),  and  8th  order  (triangle  symbols)  expansion  of  the 
atom-atom  potential.  What  is  obvious  from  the  figure  is  that  for 
oxygen-broadening  of  H2O,  calculations  utilizing  only  an  electrostatic 
expansion  of  the  intermolecular  potential  are  inadequate.  There  is 
an  improvement  in  the  agreement  between  experiment  and  theory 
when  the  4th  order  expansion  is  used.  Finally,  when  calculations  are 
made  with  the  8th  order  expansion  of  the  atom-atom  potential  the 
best  agreement  with  measurement  is  obtained.  We  comment  that 
limitations  in  our  current  computational  facilities  limits  us  to  the  8th 
order  expansion  for  this  system.  Similar  results  are  observed  for 
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calculations  of  the  line  shift  for  these  transitions.  In  Figure  11,  the 
measured  line  shifts  from  Grossmann  and  Browell119  with  error  bars 
and  the  0th,  4th,  and  8th  order  expansion  of  the  atom-atom  potential 
are  plotted  versus  the  energy  ordered  index  with  the  same  symbols 
as  Figure  10.  There  is  poor  agreement  between  measurement  and 
calculation  with  the  0th  order  calculation,  going  to  best  agreement  for 
the  8th  order  calculation. 


Figure  10  Halfwidth  in  units  of  cnr1  atm'1  for  water  vapor 
transitions  in  the  3vi+V3  vibrational  band  perturbed  by  O2 
versus  energy  ordered  index  (see  text).  Measurements  of 
Grossmann  and  Browell119  are  plotted  with  error  bars, 
0th  order  calculations  with  diamond  symbols,  4th  order 
calculations  with  circle  symbols,  and  81*1  order 
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calculations  with  triangle  symbols  (order  refers  to  the 
expansion  of  the  atom-atom  potential). 


Figure  11  Line  shifts  in  units  of  cm*1  atm'1  for  water  vapor 
transitions  in  the  3vi+V3  vibrational  band  perturbed  by  O2 
versus  energy  ordered  index  (see  text).  Measurements  of 
Grossmann  and  Browell119  are  plotted  with  error  bars, 
Oth  order  calculations  with  diamond  symbols,  4th  order 
calculations  with  circle  symbols,  and  8th  order 
calculations  with  triangle  symbols  (order  refers  to  the 
expansion  of  the  atom-atom  potential). 

The  importance  of  the  atom-atom  potential  varies  according  to 
the  strength  of  the  electrostatic  interaction.  In  general,  we  have 
found  that  for  systems  for  which  the  moments  of  the  molecules  are 
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large,  e.g.  CO2,  the  atom-atom  potential  has  a  very  small  effect  on  the 
halfwidth  and  line  shift  (<1%).  For  nitrogen  as  a  perturbing  species, 
the  effect  is  small  for  many  of  the  transitions,  however,  there  are 
some  transitions  where  the  effect  can  be  15%  of  the  calculated 
halfwidth.1 17  It  is  found  that  the  transitions  where  the  atom-atom 
potential  has  an  effect  can  be  classified  as  weak-collision  systems  in 
the  definition  of  Oka.105  Similar  results  are  found  for  hydrogen  as 
the  perturbing  species.160  For  O2  as  the  perturbing  species,  the  high- 
order  expansion  is  necessary  to  get  any  reasonable  agreement  with 
experimental  values.  In  the  calculations  presented  here,  we  have 
expanded  the  atom-atom  potential  to  8l1*,  8th,  4^,  and  4th  order  for 
N2,  O2,  CO2,  and  H2  as  the  perturbing  species,  respectively. 

Next,  we  focus  on  the  line  shifts.  It  is  observed  that  the  line 
shifts  for  water  vapor  are  both  positive  and  negative  depending  on 
the  vibrational  band  in  question.  This  can  be  understood  in  terms  of 
the  components,  Si  and  /m{S2}  that  contribute  to  the  line  shift.116 
The  Sj  term  depends  on  the  vibrational  dependence  of  the 
polarizability,  Eq.  (36).  Utilizing  the  coefficients  of  Ref.  140  this 
contribution  is  negative  and  will  increase  with  the  number  of  quanta 
in  the  vibrational  state.  The  imaginary  part  of  the  S2  term  can  be 
either  positive  or  negative.  For  transitions  for  which  the  change  in 
the  vibrational  quantum  numbers  is  small,  the  lines  shifts  will  be 
both  positive  and  negative.  As  the  change  in  the  number  of 
vibrational  quanta  for  the  transition  increases,  the  line  shifts  become 
increasingly  more  negative.  Thus,  for  the  transitions  considered 
here,  the  line  shifts  for  the  3V1+V3  and  the  2V]  +  2V2  +  V3  bands  are 
negative  and  those  for  the  Vi,  2V2,  and  V3  bands  are  both  positive 
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and  negative.  This  is  demonstrated  in  Figures  12  and  13  for  N2- 
induced  pressure  shifts  in  the  3V1  +  V3  band  of  H2O  and  for  CO2- 
induced  pressure  shifts  in  the  Vi,  2V2,  and  V3  bands  of  H2O  , 
respectively.  In  Figure  12,  the  N2  induced  pressure  shifts  in  the 
3V1+V3  band  of  H2O  are  plotted  versus  the  energy  ordered  index. 
Shown  are  the  experimental  data  of  Grossmann  and  Browell119  with 
error  bars,  the  RS2  calculated  values  are  the  diamond  symbols,  the 
RSi  calculated  values  are  the  circle  symbols,  and  the  values  from  the 
CRBF  calculation  are  triangles.  What  is  apparent  is  that  as  the 
imaginary  terms  are  included,  agreement  with  experiment  improves 
dramatically.  Note,  the  line  shifts  from  the  RS2  calculation  are  both 
positive  and  negative  and  the  line  shifts  from  the  RS]  calculations  are 
all  negative,  as  discussed  above.  In  Figure  13,  the  measured  line 
shifts  of  Gamache  et  al.116  with  error  bars  are  plotted  versus  the 
transition  number  (see  Table  I  of  Ref.  116).  Included  in  the  plot  are 
the  CRBF  calculated  values  as  solid  circles,  the  RS2  calculated  values 
as  solid  triangles,  and  the  RSi  values  as  solid  squares.  For  transitions 
numbered  17  and  27  there  are  no  measured  values  thus  only  the 
calculated  values  are  plotted.  The  RSi  values  are  all  small  due  to  the 
number  of  vibrational  quanta  involved  in  these  transitions  and  all 
negative  as  explained  above.  Thus,  for  this  system,  the  line  shift  is 
mostly  due  to  the  imaginary  part  of  S2  with  small  corrections  coming 
from  the  Si  terms  resulting  in  line  shifts  that  are  both  positive  and 
negative.  The  agreement  between  the  experimental  results  and  the 
CRBF  calculated  values  is  good  with  most  of  the  calculated  values 
falling  within  the  error  estimate  of  the  measurements.  This  is 
further  summarized  below. 


Line  Shift  / 


Figure  12  Line  shifts  in  units  of  cm1  atm-1  for  water  vapor 
transitions  in  the  3vj+V3  vibrational  band  perturbed  by 
nitrogen  versus  energy  ordered  index  (see  text). 
Calculations  with  various  levels  of  the  imaginary  terms: 
RS2  calculated  values  (diamond  symbols),  RS]  calculated 
values  (circle  symbols),  CRBF  calculated  values  (triangle 
symbols)  compared  with  the  experimental  data  of 
Grossmann  and  Browell.119 
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Figure  13  Line  shifts  for  H2O  transitions  broadened  by  CO2. 

Plotted  are  line  shifts  in  units  of  cm*1  /atm  versus 
transition  number  from  Table  1  of  Ref.  116; 
Measurements  with  error  bars,  CRBF  as  solid  circles, 
RS2  as  solid  triangles,  and  RSi  as  solid  squares. 


Last,  we  focus  on  the  imaginary  contributions  to  the  halfwidth. 
As  noted,  the  complex  Robert-Bonamy  formalism  leads  to  the 
imaginary  components  of  the  expanded  S  matrix  having  a 
contribution  to  the  halfwidth.  The  contribution  of  these  components 
to  the  halfwidth  has  remained  a  relatively  unexplored  facet  of  the 
theory  because  some  of  the  imaginary  resonance  function  integrals 
have  only  been  derived  recently.115  Previous  analysis  of  these 
contributions  for  nitrogen  and  oxygen  broadening  of  water  vapor 
transitions114’115’117’118  indicated  that  the  imaginary  parts  do 
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contribute  to  the  halfwidth.  For  the  H2O-N2  system,  the  difference 
between  the  RRB  and  CRBF  calculations  are  generally  small  with 
exceptions  up  to  10-15%  observed.  For  the  H2O-O2  system,  the 
importance  of  the  imaginary  terms  is  more  pronounced,  accounting 
for  up  to  25%  difference  between  RRB  and  CRBF  calculations.  This  is 
shown  in  Figure  14  where  the  halfwidth  is  plotted  versus  the  energy 
ordered  index  for  transitions  in  the  3V1+V3  band  of  H2O  with  O2  as 
the  perturbing  gas.  Plotted  are  the  measurements  of  Grossmann  and 
Browell119  with  error  bars,  the  RRB  calculated  values  as  circle 
symbols,  and  the  CRBF  calculated  results  as  triangle  symbols.  In 
almost  all  cases,  there  is  a  large  difference  between  the  RRB  and 
CRBF  calculated  values  with  the  CRBF  values  demonstrating  good 
agreement  with  measurement.  Analysis  indicates  that  the 
transitions  for  which  the  imaginary  terms  make  a  contribution  to  the 
halfwidth  are  also  those  for  which  the  atom-atom  terms  are 
important  (intermediate  to  weak  collisions105). 

The  H2O-CO2  system  is  characterized  by  strong  collisions,  thus 
one  might  assume  that  the  imaginary  terms  are  unimportant  in  the 
calculation  of  the  halfwidth.  This,  however,  is  not  what  is  observed 
from  the  calculations.116  In  Figure  15  the  measured  halfwidths  with 
error  bars  from  Ref.  116  are  plotted  versus  the  transition  number 
from  Table  I  of  Ref.  116  with  the  calculated  CRBF  values  as  solid 
circles,  the  RS2  values  as  solid  triangles,  and  the  RRB  values  as  solid 
squares.  From  the  figure,  it  is  apparent  that  as  the  imaginary 
components  are  included  in  the  calculation  better  agreement  with 
experiment  is  obtained.  The  difference  between  the  RRB  and  CRBF 
calculation  of  the  halfwidth  is  greater  than  the  desired  accuracy  of 


the  spectroscopic  community89,90  for  more  that  two-thirds  of  the 
lines  studied116  This  is  demonstrated  in  Figure  16  where  the  percent 
differences  between  the  RRB  and  CRBF  calculations,  (CRBF-RRB)/CRBF, 
are  shown  versus  the  transition  number  (Table  I  of  Ref.  116).  For 
hydrogen  broadening  of  H2O,  we  have  observed  that  the  imaginary 
contribution  to  the  halfwidth  is  very  small.160 


Figure  14.  Halfwidths  in  units  of  cnr1  atm-1  for  water  vapor 
transitions  in  the  3vj  +  V3  vibrational  band  perturbed  by 
oxygen.  Measurements  of  Grossmann  and  Browell119  are 
plotted  with  error  bars,  results  by  RRB  calculations  with 
circle  symbols,  and  results  by  CRBF  calculations  with 
triangle  symbols. 


Half  width  /  cm"1  atm" 


Transition 

Figure  1 5.  Halfwidths  for  H2O  transitions  broadened  by  CO2.  Plotted 
are  halfwidths  in  units  of  cm'1  atm"!  vs.  transition 
number  from  Table  I  of  Ref.  116;  Measurements  with 
error  bars,  CRBF  as  solid  circles,  RS2  as  solid  triangles,  and 
RRB  as  solid  squares. 
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Figure  16.  Percent  difference  between  CRBF  and  RRB  calculations, 
(CRBF-RRB)/CRBF,  vs.  transition  number  from  Table  I  of 
Ref  11 6. 


To  summarize,  we  have  made  CRBF  calculations  of  the 
pressure-broadened  halfwidth  for  772  transitions  of  H2O  perturbed 
by  N2, 02,  CO2,  and  H2,  and  compare  these  results  with  measurement 
These  data  are  for  the  rotational,  vi,  V2,  V3,  2v2,  3vi+V3  and  2vi+2v2+v3 
bands  of  H2O.  In  addition,  pressure-induced  line  shifts  for  212 
transitions  of  water  vapor  perturbed  by  N2,  O2,  and  CO2  are 
calculated  by  the  CRBF  method  and  compared  to  measurements. 
Tables  16-19  report  the  average  percent  differences  between  the 
CRBF  calculated  halfwidth  and  line  shift  and  the  measured  values 
from  Refs.  119,  119,  116,  and  122,  respectively.  For  hydrogen  as  the 
perturbing  species  there  are  enough  data  to  calculate  a  meaningful 


77 


standard  deviation  and  this  is  included  in  Table  19.  We  have 
calculated  the  average  percent  difference  between  the  measured  and 
calculated  line  shifts  for  the  3vi+v3  and  2vi+2v2+v3  bands  of  H2O  but 
not  for  the  vi,  2v2,  and  V3  bands.  This  is  because  in  the  former  case 
the  line  shifts  are  large  enough  in  magnitude  for  a  percent  difference 
to  be  useful.  For  the  later  bands,  the  line  shifts  are  relatively  small 
in  magnitude,  thus,  we  are  making  comparisons  of  small  numbers. 
In  fact,  for  these  data  the  experimental  uncertainty,  as  small  as 
0.002  cm*1 /atm,  can  be  as  large  as  100%  of  the  line  shift.  Thus,  it 
does  not  make  sense  to  talk  of  a  percent  difference  for  these 
comparisons.  For  the  pressure-broadened  halfwidths,  the  overall 
average  percent  difference  for  the  772  comparisons  is  -0.8  percent. 
The  average  percent  difference  for  the  183  comparisons  made  for 
the  line  shift  in  Tables  16  and  17  is  -3.6%.  Perhaps  it  is  better  to 
note  that  for  the  line  shifts,  the  CRBF  calculated  values  generally  fall 
within  or  are  close  to  the  experimental  uncertainties. 

Table  16.  Average  percent  differences  between  the  measured119 
and  CRBF  calculated  halfwidths  and  line  shifts  for  the 
water  vapor-nitrogen  system. 


Band 

Parameter 

#  lines 

APDa 

AAPDb 

3vi+  V3 

7 

49 

-2.3 

3.4 

2vi+  2v2+  V3 

y 

22 

-7.5 

10.9 

3vi+  V3 

5 

55 

-4.1 

1 1.0 

2vj+  2v2+  V3 

8 

25 

-13.3 

18.0 

a  Average  percent  difference 
b  Average  absolute  percent  difference 
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Table  17. 


Average  percent  differences  between  the  measured11 9 
and  CRBF  calculated  halfwidths  and  line  shifts  for  the 
water  vapor-oxygen  system. 


Band 

Parameter 

#  lines 

APDa 

AAPDb 

3vi+  V3 

Y 

76 

-0.2 

4.4 

2vi+  2v2+  V3 

Y 

28 

-0.6 

4.9 

3vi+  V3 

8 

78 

7.8 

8.4 

2vi+  2v2+  V3 

8 

25 

-28.2 

28.2 

a  Average  percent  difference 
b  Average  absolute  percent  difference 


Table  18.  Average  percent  differences  between  the  measured116 
and  CRBF  calculated  halfwidths  and  line  shifts  for  the 
water  vapor-carbon  dioxide  system. 


Band 

Parameter 

#  lines 

APDa 

vi 

Y 

8 

0.4 

2v2 

Y 

1 

-0.8 

V3 

Y 

22 

1.8 

a  Average  percent  difference 


Table  19. 


Average  percent  differences  between  the  measured122 
and  CRBF  calculated  halfwidths  and  line  shifts  for  the 
water  vapor-hydrogen  system. 


Band 

Parameter 

#  lines 

APDa 

SDb 

rotation 

Y 

64 

-1.9 

4.2 

vi 

Y 

92 

2.7 

4.6 

V2 

Y 

273 

-2.8 

5.4 

V3 

Y 

137 

2.1 

4.2 

a  Average  percent  difference 
b  Standard  deviation 


In  conclusion,  we  give  the  following  comments  with  respect  to 
calculations  of  pressure-broadened  halfwidths  and  pressure-induced 
line  shifts  for  ro-vibrational  transitions  of  water  vapor.  The  effects 
of  the  imaginary  components,  the  choice  of  potential,  and  the  order 
of  expansion  of  the  potential  vary  from  transition  to  transition.  To 
insure  a  reasonably  correct  calculation,  the  CRBF  approach  with  the 
leading  terms  of  the  electrostatic  potential,  the  atom-atom  potential 
expanded  to  at  least  4th-order,  and  the  isotropic  induction  and 
dispersion  potentials  should  be  used. 

The  results  of  these  calculations  demonstrate  a  profound 
dependence  of  the  halfwidth  on  the  imaginary  terms  in  the 
calculation  for  some  of  the  perturbing  species  considered  here.  For 
the  772  transitions  studied  for  4  perturbing  species,  the  CRBF 
calculation  of  the  halfwidth  is  in  better  agreement  with  experiment 
than  calculated  values  based  on  other  theoretical  methods.  Almost 
always,  as  the  imaginary  components  are  added  and  the  atom-atom 
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potential  is  included  to  high  order,  the  results  move  in  the  direction 
of  measurement.  And  last,  the  line  shifts  calculated  by  the  CRBF 
method  agree  very  well  with  experiment. 

The  CRBF  formulation  offers  a  cutoff-free,  scaling-parameter- 
free  approach  to  calculating  the  halfwidth  and  line  shift.  The  authors 
are  testing  the  CRBF  formulation  on  other  radiator  and  perturbing 
species  and  investigating  the  temperature  dependence  of  the 
halfwidth  and  line  shift.  This  work  will  be  presented  in  a  future 
publication.  It  appears  that  the  CRBF  formalism  is  close  to  satisfying 
the  spectroscopic  community’s  need  to  know  such  parameters  to  an 
accuracy  of  1  -5%. 89. 90 

4.0  PROGRAMS  FOR  THE  HITRAN  DATABASE 

A  new  version  of  the  BANDSUM  program  was  written  which 
works  with  the  new  .CMN  files  etc.  Also  BD-VIBS.FOR  was  updated 
again  for  the  Schuman-Runge  band  information.  We  have  added  v=2 
of  X3E'g  and  v=0  to  19  of  B3£-g.  The  list  of  vibrational  levels  of 
Brown  and  Urban  have  been  added  to  BD-VIBS.FOR. 

A  number  of  supplemental  programs  that  will  accompany  the 
database  have  been  written.  These  are  programs  that  calculate  the 
natural  isotopic  abundance,  the  molecular  masses,  programs  to  check 
the  band  center  list,  and  similar  types  of  programs.  All  programs 
have  been  delivered  to  Dr.  L.  Rothman  for  use  in  the  maintenance  of 
the  database.  In  addition,  updates  to  several  of  the  BLOCK  DATA 
routines  were  made. 


A  new  table  has  also  been  added  to  the  database  which 
contains  the  information  for  the  isotopic  label,  natural  isotopic 
abundance,  the  partition  sum  at  296K,  the  state-independent  nuclear 
degeneracy  factor,  and  the  molar  mass  for  each  species.  This  table  is 
called  MOLPARAM.TXT  and  is  distributed  with  the  database. 

4.1  Transition  Moment  Squared 

The  calculation  of  the  transition  moment  squared  for  all  species 
on  the  database  was  reviewed.  Since  the  last  version  of  the  database 
there  have  been  many  changes  to  the  vector  indices  that  control  the 
database.  These  changes  had  not  been  made  to  the  TR-MOM 
program.  In  addition  to  these  changes,  there  were  many  cases  where 
the  transition  moment  squared  could  not  be  calculated  due  to  the 
partition  sums  not  being  present.  For  the  1995  database,  there  is 
also  the  addition  of  several  new  molecules  and  additional 
isotopomers  for  some  molecules.  All  of  these  factors  need  to  be 
included  in  the  1995  version  of  TR-MOM. 

The  first  changes  were  to  add  to  the  code  the  correct  .CMN  files 
in  order  to  have  correspondence  with  the  latest  molecule,  abundance, 
isotopomer  and  vibrational  labeling.  Next,  subroutines  were  written 
for  all  the  new  molecular  species  that  have  been  added  to  the 
database.  The  old  subroutines  were  updated  for  new  isotopomers 
where  appropriate.  The  subroutines  that  defaulted  a  value  of  R^ 
equal  to  zero  were  also  corrected.  The  individual  subroutines  were 
checked  for  consistent  use  of  partition  sums,  degeneracy  factors,  etc. 
The  updates  to  the  transition  moment  squared  program  that  were 
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discussed  in  the  last  quarterly  report  were  implemented.  These 
include  the  vector  indices  (molecule  and  isotope  codes),  isotopic 
abundance,  updated  partition  sums  and  degeneracies,  and  other 
corrections  in  the  codes.  One  major  improvement  is  the  ability  to 
calculate  R2  for  all  the  species  on  the  main  database.  Thus,  the 
subroutines  that  defaulted  to  a  zero  value  of  R2  now  correctly 
calculate  the  transition  moment.  For  the  molecular  species  new  to 
the  database  subroutines  were  written  to  calculate  R2.  For  new 
isotopomers  old  subroutines  were  updated  to  consider  the  new 
species.  The  individual  subroutines  were  checked  for  consistent  use 
of  partition  sums,  degeneracy  factors,  etc.  Test  calculations  unveiled 
some  additional  errors  which  were  corrected  and  the  new  version  of 
the  TR-MOM  program  was  linked  with  the  latest  versions  of  BD-ISO, 
BD-MOL,  BD-VIBS,  BD-ABUN,  and  BD-QT. 

4.2  State  Independent  Degeneracy  Factors 

Table  20  lists  the  molecule  number,  molecule,  isotopomer  code, 
species  code,  and  the  state  independent  degeneracy  factors  that  are 
used  for  determining  the  partition  function,  correct  transition 
moment  squared  factor,  and  conversion  between  intensity,  R2, 
Einstein-A  coefficients,  etc. 
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Table  20 


State  Independent  Degeneracy  Factors 


1 

h2o 

161 

1 

1 

181 

2 

1 

12  HNO3 

146 

45 

6 

171 

3 

6 

162 

4 

6 

13  OH 

61 

46 

r 

81 

47 

r 

2 

co2 

626 

5 

1 

62 

48 

r 

636 

6 

2 

628 

7 

1 

14  HF 

19 

49 

4 

627 

8 

6 

638 

9 

2 

15  HC1 

15 

50 

8 

637 

10 

12 

17 

51 

8 

828 

11 

1 

728 

12 

6 

16  HBr 

19 

52 

8 

11 

53 

8 

3 

03 

666 

13 

1 

668 

14 

1 

17  HI 

17 

54 

12 

686 

15 

1 

668 

16 

1 

18  CIO 

56 

55 

r 

686 

17 

1 

76 

56 

r 

667 

18 

6 

676 

19 

6 

19  OCS 

622 

57 

1 

624 

58 

1 

4 

n2o 

446 

20 

9 

632 

59 

2 

456 

21 

6 

822 

60 

1 

546 

22 

6 

623 

61 

4 

448 

23 

9 

447 

24 

54 

20  H2CO 

126 

62 

1 

136 

63 

.  2 

5 

CO 

26 

25 

1 

128 

64 

1 

36 

26 

2 

28 

27 

1 

21  HOC1 

165 

65 

8 

27 

28 

6 

167 

66 

8 

38 

29 

2 

37 

30 

12 

22  N2 

44 

67 

1 

6 

CH4 

211 

31 

1 

23  HCN 

124 

68 

6 

311 

32 

2 

134 

69 

12 

212 

33 

3 

125 

70 

4 

7 

o2 

66 

34 

1 

24  CH3CI 

215 

71 

4 

68 

35 

1 

217 

72 

4 

67 

36 

6 

25  H202 

1661 

73 

4 

8 

NO 

46 

37 

r 

56 

38 

r 

26  C2H2 

1221 

74 

1 

48 

39 

r 

1231 

75 

8 

9 

so2 

626 

40 

1 

27  C2H6 

1221 

76 

64 

646 

41 

1 

1111 

77 

2 

28  PH3 

10 

no2 

646 

42 

3 

29  COF2 

269 

78 

1 

11 

nh3 

4111 

43 

3 

5111 

44 

2 
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Table  20. 

State 

Independent 

Degeneracy 

Factors 

30  SF6 

29 

79 

1 

7646 

87 

12 

31  H2S 

121 

80 

1 

36  NO+ 

46 

88 

3 

141 

81 

1 

131 

82 

4 

37  HOBr 

169 

89 

8 

161 

90 

8 

32  HCOOH 

126 

83 

4 

38  H4C2 

1221 

91 

1 

33  HO2 

166 

84 

2 

1231 

92 

32 

34  O 

6 

85 

NA 

35  CIONO2 

5646 

86 

12 

5.0  SPECTRAL  PARAMETERS  FOR  SPECIES  ON  THE  DATABASE 

5.1  Errors  in  the  1992  data  for  Hydrogen  Halides  and  Carbon 
Monoxide 

Dr.  Curtis  Rinsland,  of  NASA  Langley  Research  Center,  reported 
to  us  that  there  were  problems  with  the  HC1  fundamental  (l*-0)  line 
positions  compared  with  measurements  he  has  made.161  While 
investigating  this  problem,  it  was  suspected  that  there  were  errors  in 
the  other  hydrogen  halides  and  CO.  For  HC1,  comparisons  were  made 
with  the  data  of  Rinsland,  the  1986  HITRAN  database  data162  and 
the  1992  data.163  This  study  revealed  errors  in  the  data  for  the  HC1 
fundamental  band  on  the  1992  database.3  Further  investigation 
showed  the  data  for  all  bands  of  HC1  on  HITRAN92  have  some  errors. 

These  problems  led  us  to  question  whether  the  data  for  the 
other  closed  shell  diatomics  are  correct.  The  reference  point  was 
taken  as  the  1986  HITRAN  database.  The  line  positions  (HITRAN92- 
HITRAN86)  and  the  ratio  in  the  line  intensities 
(HITRAN92/HITRAN86)  were  compared  for  HF,  HCI,  HBr,  HI,  and  CO. 
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The  comparisons  were  made  and,  with  the  exception  of  HC1 

positions,  the  results  indicated  that  the  data  for  the  other  diatomics, 
positions  and  intensities  (HC1  included),  are  not  grossly  different, 
except  for  several  typographical  errors.  These  errors  were 
summarized  and  delivered  to  Dr.  L.  Rothman. 

For  the  line  positions,  we  have  calculated  the  energies  of  the 
diatomic  molecules  using  the  Dunham  expansion  and  compared  the 

values  with  the  1986  and  1992  databases.  This  indicated  that  the 

1992  data  was  in  error.  From  discussions  with  the  author  of  this 

data,163  we  discovered  the  data  files  were  typed  by  a  secretary  from 
computer  output  and  concluded  this  to  be  the  source  of  the 

typographical  errors.  This  data  needs  to  be  corrected  for  the  next 
version  of  the  database. 

5.2  Vibrational  Band  Centers 

We  have  modified  the  file  of  vibrational  band  centers  to  be 

more  easily  maintainable.  The  existing  file  had  information  about 

the  molecule,  isotopic  species,  vibrational  band,  band  center,  and  a 
written  reference  in  an  A40  format.  It  is  our  intentions  to  update 
this  file  with  improved  band  centers.  What  we  have  done  is  to 

reformat  the  file  so  that  it  can  be  read  in  free  format  and  to 

eliminate  the  need  for  writing  the  reference.  The  reference  is  now 

the  HITRAN  reference  code  and  can  be  found  in  TABLES. 3  of  the 

database.  We  have  gone  through  the  file  and  cross-referenced  with 
TABLES. 3.  For  the  references  that  existed  in  TABLES. 3  we  must  use 
the  same  numbers,  for  others  we  added  numbers  and  references  to 
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TABLES. 3  accordingly.  The  new  file  of  band  centers  has  the  structure 
Mol  number,  HITRAN  isotope  code,  spectroscopic  vibrational  quanta 
(’  <-  "),  the  band  center  in  cm'1,  reference  number. 

5.3  Other  Species  Data  for  HITRAN 

New  molecules  and  isotopic  species  have  been  added  to  the 
database.  The  number  of  molecules  is  now  36  (  up  from  32  on  the 
1992  Edition).  The  additions  include  the  667  and  676  isotopic 
species  of  ozone  and  the  37  isotopomer  of  CO.  The  new  additions 
with  the  corresponding  molecule  number  a  given  in  Table  21. 

The  hydroperoxyl  radical,  HO2  has  been  added  as  molecule  33, 
it  is  an  asymmetric  rotor  and  is  vibrational  Class  6  and  in  rotation 
Group  1.  The  oxygen  atom  is  molecule  34,  it  has  no  vibrational 
structure  (an  atom)  and  is  a  special  case.  It  has  been  put  in 
vibrational  Class  1  and  rotational  Group  2.  Next,  molecule  35,  is 
CIONO2,  this  has  been  placed  with  the  other  large  species,  i.e. 
vibrational  class  10  and  in  rotational  Group  1.  NO+  has  been  added 
as  molecule  #  36  and  is  in  vibrational  Class  1. 

The  new  species  have  the  following  isotopic  species  codes;  one 
isotope  of  HO2  is  present  labeled  166  for  H16C>2,  oxygen  atom  (160) 
has  the  code  6,  for  CIONO2  there  are  two  isotopes  labeled  5646  for 
35C116014N1602  and  7646  for  37cil6014N1602. 

There  have  also  been  additions  to  the  vibrational  tables  for 
some  molecules.  In  particular,  the  O2  table  has  been  extended  to  29 
elements  with  the  additions  of  the  vibrational  state  labels  "  BO"  to 
"  B19".  Dr.  Linda  Brown  has  extended  the  list  for  NH3,  this  still 
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needs  to  be  reviewed.  A  complete  list  can  be  found  in  the  routine 
BD-VIBS.FOR. 


Table  21  New  Molecules  and  Isotopomers  Added  to  the  HITRAN 
Database  since  the  1986  Edition. 


Molecule 

Molecule  Code 

Isotope  Code 

16016C)170 

3 

667 

160170160 

3 

676 

16Q12033S 

1  9 

623 

cof2 

29 

269 

sf6 

30 

29 

h2s 

3  1 

121 

3  1 

141 

3  1 

1  3  1 

HCOOH 

32 

126 

HO 

33 

166 

0 

34 

6 

CIONO 

35 

5646 

7646 

N0+ 

36 

46 

HOBr 

37 

169 

161 

h4c2 

38 

1221 

1231 

With  the  new  additions,  many  changes  had  to  be  made  to  the 
block  data  routines  that  accompany  most  of  the  HITRAN  programs. 
A  new  structure  was  adopted  for  these  routines.  We  now  use 
INCLUDE  statements  for  abundance,  isotopomers,  Q(T),  species,  and 
VIBS.  All  routines  were  updated  with  the  new  information. 


5.4  Ozone  Data  for  the  1996  Edition  of  HITRAN 


Nine  files  of  ozone  data  for  the  1994  Edition  of  HITRAN  were 
obtained  from  Dr.  C.  P.  Rinsland.164  The  data  contained  in  these  files 
are  listed  in  Table  22.  The  files  were  decompressed  and  run  through 
the  ozone  halfwidth  program  which  adds  the  air-  and  self-broadened 
halfwidths  and  the  air-broadening  temperature  exponent  to  the  files 
along  with  the  reference  codes.  The  data  was  then  transferred  to  Dr. 
Rothman  for  addition  to  the  1996  HITRAN  database. 

5.5  NH3  for  HITRAN96 

New  and  updated  data  were  obtained  for  NH3  from  L.  Brown 
and  S.  Urban.  The  changes  in  the  data  include  new  calculated 
wavenumber,  intensity,  lower  state  energies,  and  self-  and  air- 
broadened  halfwidths,  the  notation  for  the  rotational  quantum 
numbers  has  been  changed,  and  many  new  vibrational  bands  have 
been  included.  The  data  sent  by  both  Urban  and  Brown  has  used  a 
vibrational  labeling  scheme  that  does  not  agree  with  the  HITRAN 
system  and  must  be  changed  to  be  compatible  with  the  database. 
The  operations  were  made  on  the  original  five  data  files.  These  are 
listed  in  Table  23.  The  HINH3.JPL  file  was  data  of  Brown. 
NH3_3000.DAT  and  NH3_4000.DAT  were  the  new  data  of  Urban,  and 
that  NH3_HITR.DAT  was  an  incomplete  transfer  of  the  file 
NH3_2500.DAT.  Furthermore,  the  file  NH3_2500.DAT  is  the  1992 
HITRAN  data  for  NH3  in  the  region  0-2154  cm'1.  This  file  was 
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generated  by  Urban  and  contains  the  new  vibrational  assignment 
codes  and  updated  rotational  quantum  number  notation.  Other  than 
the  quantum  codes  the  data  match  HITRAN92. 


Table  22.  Ozone  Files  for  the  1996  HITRAN  Compilation 


File 

ISO 

V’ 

V" 

Status 

Region 

#  lines 

HIT03941.ARC 

666 

[(004),  (103), 

(000) 

Nf 

2.5pm 

4050 

(310)] 

HIT03942.ARC 

666 

[(013),  (112)] 

(000) 

N 

2.7pm 

2221 

HIT03943.ARC 

666 

[(004),  (103), 

(010) 

N 

3.3pm 

476 

(310)] 

HIT03944.ARC 

666 

|(003),  (102), 

(000) 

Rtt 

3.3pm 

9746 

(201),  (300)] 

HIT03945.ARC 

666 

[(003),  (102), 

[(100), 

R 

4.8pm 

11262 

(201),  (300)] 

(001)] 

HIT03946.ARC 

666 

[(003),  (102), 

[(200), 

N 

10.0pm 

4287 

(201),  (300)]  (101),  (002)] 

HIT03947.ARC 

666 

|(004),  (103), 

[(100), 

N 

3.3pm 

5179 

(310)] 

(001)] 

1(013), (112)] 

(010) 

HIT03948.ARC 

667 

(010) 

(000) 

N 

14pm 

5641 

HIT03949.ARC 

676 

(QiQ) _ 

(000) 

N 

14pm 

2846 

t  New  data,  not  present  on  HITRAN  1992. 
tt  Replace  data  that  are  on  HITRAN  1992. 


Table  23 


NH3  Data  Files 


File 

#  Lines 

Wmin^ 

NH3_3000.DAT 

836 

3118 

3675 

NH3_2500.DAT 

4855 

0 

2154 

NH3_4000.DAT 

1130 

4036 

4710 

NH3_H1TR.DAT 

2570 

0 

1  249 

HINH3.JPL 

3370 

2150 

5295 

t  in  cnr1 


The  data  was  grouped  into  2  files,  BROWN. PA  1  and  URBAN.PA1. 
The  first  is  the  file  H1NH3.JPL  and  the  second  is  NH3_2500.DAT, 
NH3_3000.DAT,  and  NH3_4000.DAT.  A  program  was  written  to  read 
these  data  and  make  the  necessary  corrections  to  the  vibrational 
indices,  add  references  and  error  codes.  The  data,  which  were  sent 
by  L.  Brown  to  L.  Rothman,  were  to  have  updated  halfwidths  already. 
Upon  inspection,  it  was  found  that  the  JPL  file  had  incorrect 
halfwidths  (changing  by  2  orders  of  magnitude  for  similar  J  K 
transitions).  Closer  inspection  revealed  that  the  Urban  data,  which 
looked  correct,  were  also  in  error  by  a  few  percent.  This  was  traced 
to  a  units  conversion  error  in  L.  Brown's  program.165  We  took  the 
data  of  A.  Pine166  for  N2-  and  02-broadening,  which  are  reported  in 
cm_1/MPa  and  converted  these  to  cnr'/atm.  The  data  were  then 
fitted  by  the  polynomial  expression  used  by  Brown1 9 

Y(J",K")  =  a()  +  arJ”  +  a2K”  +  a3J"2  +a4J"K"  .  (40) 

The  resulting  coefficients  are  given  in  Table  24.  These  are  then  used 
with  the  self-broadening  coefficients  of  Brown165  to  add  these  to  the 


two  files.  The  resulting  files  have  been  called  BROWN. PA2  and 
URBAN.PA2.  These  need  to  be  checked  and  L.  Brown  suggest  the 
V2-V2  data  on  the  database  is  too  intense  by  a  factor  of  2.  Using  a 
program  that  identifies  the  lines  of  the  V2-V2  band  (labeled  6  and  10, 
or  10  and  6)  the  intensities  for  these  lines  was  divided  by  2.  There 
were  199  lines  for  the  OlOOs-OlOOa  band  and  394  lines  for  the 
OlOOa-OlOOs  band.  After  this  the  halfwidth  program  was  run  to  add 
the  self-  and  air-broadened  halfwidths  to  the  data.  Band  sums  were 
made  on  the  final  files,  called  BROWN. PA2  and  URBAN. PA3.  It  was 
discovered  that  the  labels  used  for  certain  vibrational  states  were 
incorrect  and  inverted  in  cases.  Comparing  with  the  "bean  counting" 
program  of  L.  Brown  other  problems  arose  in  that  the  states  0000s 
and  0000a  are  both  labeled  1  and  2. 

From  Herzberg167  we  find  that  Vi  (1000-0000)  is  a  parallel 
band  as  well  as  11000-0000.  For  both  of  these  vibrational 
transitions  the  selection  rules  a-s,  and  s-a  apply.  This  indicated  that 
the  band  sums  sent  by  L.  Brown  had  the  data  labeled  incorrectly. 
Checking  Urban's  original  data  gives 

1000s -0000a  called  18  2  in  the  data 
lOOOa-OOOOs  called  19  1  in  the  data. 

These  assignments  agree  with  the  correspondence  list  that  has  been 
created  for  ammonia.  This  confirmed  that  the  data  was  correct  and 
the  labeling  error  was  only  in  L.  Brown's  "bean  counting"  program. 
The  final  data  files  are  BROWN. PA2  and  URBAN. PA3.  These  were 
delivered  to  L.  Rothman  for  the  1996  database. 
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Table  24.  Coefficients  of  Polynomial  in  Eq.  (40)  for  N2-  and  O2- 
Broadening  of  NH3 


ao 

ai 

32 

33 

a4 

NH3-N2 

NH3-02 

.11982 

.69178  10*1 

-.94290  10* 2 
-.44293  10* 2 

59746  10* 2 
.23350  10*2 

.63050  10*3 
.30747  10*3 

-.61553  10*3 
-.31585  10*3 

5.6  Transition  Moment  Squared  for  Ozone  Molecule 

Dr.  M.  A.  H.  Smith  of  NASA  Langley  Research  Center  sent  email 
to  Dr.  Laurence  Rothman  of  Phillips  Laboratory  (forwarded  to  me) 
stating  that  she  noticed  that  for  ozone  lines  in  the  990-1010  cm*1 
region  on  the  1996  H1TRAN  database  the  transition  moment  squared, 
IRI2,  is  smaller  than  the  1992  HITRAN  database  values,  yet  the  other 
parameters  for  these  lines  are  exactly  the  same.  To  investigate  this  I 
considered  the  calculation  of  the  transition  moment  squared  and 
what  changes  could  have  taken  place.  The  most  obvious  was  the 
improved  partition  sums  on  HITRAN96.  Running  TIPS92  and  TIPS96 
gives  the  results  in  Table  25.  This  yields  a  1996/1992  ratio  of 
Q96/Q92  =  0.99812. 
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Table  25. 


Partition  Functions  at  296K  for  Ozone  Molecule 


Q(296  K)  TIPS96 

Q(296  K)  TIPS92 

3481.86 

3488.42 

Since  the  transition  moment  squared  is  proportional  to  Q  this 
gives  a  difference  in  IRI2  values  of  IRI2(96)  =  0.998 12*IRI2(92).  I  have 
selected  some  ozone  lines  in  the  990-1010  cm-1  region  from  the 
1992  and  1996  databases.  These  values  are  presented  along  with 
scaled  1992  values  in  Table  26.  Thus,  it  appears  that  the  difference 
is  entirely  due  to  the  improvements  in  the  partition  functions. 


Table  26.  Transition  Moment  Squared  for  O3  Transitions  in  the 
990-1010  cm'1  region  from  the  1992,  scaled  1992, 
and  1996  HITRAN  Databases 


position  cm- 1 

IRI2(96) 

0.99812*IRI2(92) 

IRI2(92) 

919.716600 

1.699E-2 

1.700E-2 

1.703E-2 

921 .309800 

2.306E-5 

2.3057E-5 

2.301E-5 

923.02370 

1.705E-2 

1.7048E-2 

1.708E-2 

924.956000 

1.652E-2 

1.6519E-2 

1.655E-2 

931.375900 

8.026E-3 

8.02588E-3 

8.041E-3 

5.7  Halfwidths  for  Transitions  of  H2O2 


We  have  determined  two  linear  fits  to  halfwidth  data  for  H2O2. 
The  data  are  those  of  Devi  et  al.168  There  were  18  air-broadened 
transitions  studied  from  J"=3  to  15.  The  first  fit  is  to  all  of  the  data 
and  another  fit  to  the  data  averaged  as  a  function  of  J"  was  made. 
There  appeared  to  be  no  significant  difference  between  the  fits  so 
the  fit  to  all  the  data  was  chosen  as  the  final  one.  The  data  and  the 
fit  are  shown  in  Fig  17.  Although  there  is  much  variation  in  the  data, 
given  the  lack  of  more  data  this  fit  is  as  representative  as  can  be 
produced. 


Figure  17.  Linear  fit  in  J"  to  halfwidths  of  H2O2 
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The  fit  of  the  data  yields  the  following  formula  to  describe  the  air- 
broadened  halfwidths  of  hydrogen  peroxide; 

Y(cm_1  atm"1)  =  0.099057  +3.7703x10^  xj"  (41) 

5.8  State  Independent  Degeneracy  Factors  for  New  Species  on 
the  Database 

Several  new  species  have  been  added  to  the  database.  With 
the  additions,  appropriate  changes  need  to  be  made  to  a  number  of 
programs  in  order  to  maintain  the  indexing  structure  of  the  database. 
This  is  especially  true  when  one  of  the  added  species  comes  within 
the  current  list,  as  is  the  case  here.  The  new  species  are  HOBr,  C2H4, 
and  16012C33S.  The  last  of  these  species  is  in  the  old  list  so  care  must 
be  taken  to  arrange  the  indices  correctly.  For  HOBr  and  C2H4  there 
are  two  species  each,  H16079Br  and  H^O^^Br;  and  H212C2  and 
H12C13CH.  The  state  independent  degeneracy  factors  are  given  by 

gj  =  0(2Ii+l)  (42) 


where  the  product  is  over  all  nuclear  spins  that  do  not  couple  with 
rotation,  hence  they  are  independent  of  the  rotational  state.  The 
nuclear  spins  of  the  species  considered  here  are:  I(12C)=0,  I(16O)=0, 
I(1H)=l/2,  I(79Br)=3/2,  1(81  Br)=3/2,  I(13q=1/2.  For  I6()l2c33s  there 

is  no  inversion  or  rotational  symmetry  so  all  nuclei  are  included  in 
the  product,  gj  =  (2x0+1)  (2x0+1)  (2x-|+l)  =  4. 


96 


For  HOBr,  because  the  spins  do  not  change  for  the  two  species, 
gj  is  the  same  for  both  and  is  given  by  gj  =  (2x^+1)  (2x0+1)  (2x-|+l)  =  8  . 

For  hydrogen  peroxide  the  two  species  are  different.  The  first, 
H212C2>  has  C2 v  symmetry  so  that  the  spin  of  the  hydrogen  atoms 
couples  with  the  rotational  motion  to  give  odd  and  even  states 
statistical  factors  of  3  and  1,  respectively.  For  the  second  species, 
H12C13CH,  the  presence  of  the  carbon-13  breaks  this  symmetry  so 

that  the  state  independent  degeneracy  factor  is  given  by 
gj  =  (2x0+1)  (2xl+l)5  =  32  . 

6.0  EINSTEIN  A  COEFFICIENT,  INTEGRATED  BAND  INTENSITY, 

AND  DEGENERACY  FACTORS 

The  correct  use  of  the  nuclear  degeneracy  factors  for  relating 
the  Einstein  A  coefficient  to  the  integrated  band  intensity  has  posed 
a  problem  to  many  researchers.  The  reason  is  in  part  due  to  a  casual 
use  (or  misuse)  of  notation.  We  are  currently  involved  in  calculating 
the  line  intensities  for  the  X3Z'g  -»  a]Ag  band  of  molecular  oxygen 
from  the  Einstein  A  coefficient.  In  the  literature,  there  are  3 
determinations  of  the  integrated  band  intensity65-67  from  which  the 
Einstein  A  coefficient  is  deduced.  Unfortunately,  there  are  three 
proposed  sets  of  degeneracy  factors.  The  measurements  differ  each 
by  roughly  a  factor  of  2,  thus  the  overall  spread  is  a  factor  of  4.  Yet, 
the  derived  Einstein  A  coefficients  are  close  to  each  other  due  to  the 
use  of  incorrect  degeneracy  factors,  (see  Table  27) 
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Table  27. 


Measured  band  intensities,  derived  Einstein-A 
coefficients,  and  statistical  degeneracy  factors  for  the 
a1  Ag  (v=0)  <-  X3I'g  (v=0)  band. 


Reference 

Soo 

cm* 1  /(mol  ec*cnr2) 

Ae'AV->e"A"v” 

sec-1 

dz/du 

Badger  et  al.65 

3.66  x  lO'24 

2.58  x  10'4 

3/2 

Lin  et  al.66 

9.4  x  10-24 

1.3  x  10'4 

3/1 

Hsu  et  al.67 

2.1  x  10-24 

2.9  x  10'4 

3/1 

This  is  an  important  problem  since  the  Einstein  A  coefficients  are 
used  to  determine  line  intensities  in  other  models.  Thus,  application 
of  intensity  equations  will  only  generate  consistent  line  intensities  if 
the  statistical  degeneracy  factors  and  the  derived  Einstein-A 

coefficient  are  from  the  same  author. 2  The  line  intensities  for  the 
X3£-g  -»  a1  Ag  band  on  H1TRAN92  are  roughly  a  factor  of  2  too  small, 

this  is  due  to  incorrect  inversion  from  Badger  et  al.'s^S  A  to  S  (i.e.  the 

degeneracy  factors).  Suspicion  of  such  missing  factors  of  2,  and 
concerns  about  the  interpretation  of  upper  atmosphere  emissions, 
such  as  inferring  ozone  from  SME  (Solar  Mesophere  Explorer  via  the 
a  1  Ag  1.27pm  airglow)  were  expressed  by  Mlynczak  and  Nesbitt70 
(who  conjectured  a  significant  change  in  A  from  one  of  the  reported  S 
values  in  Table  5  (Hsu  et  al.»  and  by  Pendelton  et  al.71  Recent 
observations  in  the  mesosphere72  confirm  the  Badger  et  al.65  value  of 
the  Einstein-A  coefficient. 

The  correct  equations  are  discussed  in  Ref.  2,  however  a 

complete  derivation  was  not  possible  due  to  space  limitations.  Below 
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we  derive  the  relationships  between  the  Einstein  A  coefficient,  the 
integrated  band  intensity,  the  line  intensity  and  degeneracy  factors. 

The  discussion  is  tailored  to  the  oxygen  molecule,  however,  it 
can  easily  be  extended  to  other  systems  by  the  appropriate  labeling 
of  states.  For  O2,  which  has  2  low-lying  excited  electronic  states  and 
is  an  open  shell  electronic  system  in  the  ground  state  the  labels 
needed  are  an  electronic  quantum  number,  e,  lambda  doubling 
quantum  number.  A,  vibrational  quantum  number,  v,  rotational 
quantum  number,  J,  and  the  projection  of  J  on  the  internuclear  axis, 
M.  For  O2  the  following  short  hand  notation  for  the  Quantum 
numbers  will  be  used,  e=eAv.  The  standard  spectroscopic  notation  of 
double  primes  for  the  lower  state  and  a  single  prime  for  the  upper 
state  will  be  used  and  Penner's  labeling  scheme45  will  also  be  used. 
The  degeneracy  factors  given  by  the  product  of  the  degeneracy 
factors  for  the  quantized  motions,  dt|=dEdAdvdjdsym  with  de=(2S+l), 
dA=(2-6A,o),  dv  =  l,  dj=(2J+l),  and  dsym  is  one  for  the  heteronuclear 
species  and  the  (+)  states  of  the  homonuclear  species  and  zero  for  the 
(-)  states  of  the  homonuclear  species.  This  factor  accounts  for  the 
symmetry  restrictions  for  the  homonuclear  diatomics  (note,  this 
gives  an  overall  factor  of  1/2  in  the  partition  sum). 

We  begin  by  considering  non-degenerate  levels.  The  Einstein  A 
coefficient  in  sec1  units  is  given  in  terms  of  the  transition  moment 
squared  in  Debye  squared  units  as45 


_  64x4  IQ-36 


3  h 


(43) 
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In  the  absence  of  an  E  or  B  field,  the  states  labeled  by  M  are 
degenerate,  so  co(e"j"M")(e,J,M,)=W(e"j,,)(e’J')-  Ref-  45  tells  us  that 


AET-»eMJ 


X  Ae'j'M'— 


(44) 


where  dj'  is  the  degeneracy  of  states  summed  over,  here  2J7+1.  The 
tilde  over  the  A  implies  that  this  is  a  reduced  quantity.  Thus,  for  the 
problem  at  hand  we  get 


Aet->e"J"  =  7-  X  }°  36~  0>(E"J")(e,J’)  |R(e"J"M")(e'J'M')|2 

d J*  M",M'  3h 

=  ^4n4  IQ  36  co^e»j»)(e-j ) -L-  X  |R(e"J"m,,)(£'j,m')|2  (45) 

3  h  dj-  M»>M. 


The  line  intensity  (cnH/molec^cnr2)  in  terms  of  the  Einstein  A 
coefficients  for  the  non-degenerate  case  is  given  by45-1 


S(e“j"M“)(£'J'M',  =  f-  Ne-j-m-  — — 1 - (1-  e‘hv<£"J"M"H£'J'MVkT)  A£  j-M' 


Bn 


V(E'TM")(eTM') 


(46) 


where  n  is  the  frequency  in  sec1,  c  the  speed  of  light.  In  what 
follows  we  adopt  a  shorthand  label  for  the  radiation  field  part  of  Eq. 
(46),  i.e.  (l-  e'hv(e"J"M'')(£'J'M’)/kT^frad].  For  the  no-field  case,  the  M 

states  are  degenerate  so  what  is  measured  is  the  sum  over  these 
states. 
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S(£”j")(e'J')  - 


_  i  v  [rad] 


%KC  M",M'  ®?e"J")(e'J" 


NE"j"M"  .  (47) 


The  population  factors  can  be  written 


.(Lls**!  .(2j+I)njm  (48) 
Q(T)  Q(T) 


which  implies  that 


N£jm  =  ,  and  dejM  = 

(2J+1 )  (2J+1 ) 


JLlL 


(49) 


Using  Eq.  (49),  the  line  intensity  can  be  written 


C  _  1  Ne-r  [rad]  v  a 

^(e’T'XeT)  —  ~  2  k’T'M" 

°n c  dj-  co(e-T)(eT)  m\M’ 


(50) 


From  the  definition  of  the  reduced  Einstein  A  coefficient,  Eq.  (44)  we 
get 


S(eT)(E'J')  -  — Ne"j — — -  [rad]  Aej'^E'j"  (51) 

87CC  W(E"J”)(E'J') 

which  is  equivalent  to  Eq.  2-17  of  Penner^S  for  the  case  of  M 
degeneracy.  From  Eqs.  (45)  and  (51)  we  can  write  the  following 
result 
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S(e"J")(£T) 


=  J_  Ne-r - 1 - [rad]  ^  ^7t4  1Q:3^ 

87tc  coJE"j..)(eT)  dj"  3  h 


x  w(e”J")(e'J') 


dj' 


I 


|R 


=  —  10  36  Ne"j"  C0(e"j")(e,j')  [rad]  -J—  £  |R(e"J"M")(eTM')|2  -(52) 

3 he  dj-  M..  M. 


This  is  Eq.  (38)  of  Gamache  and  Rothman1  and  Eq.  (1)  of  Ref.  2 
without  the  isotopic  abundance  factor  used  in  HITRAN. 

We  now  wish  to  develop  the  integrated  band  intensity,  SE-,  for 
an  electronic-vibrational-rotational  band.  In  units  of 

cm*  Vmolec'cnv2  this  is  given  by 


Se,  _  871  10  36.  (0(E“)(£,)  Ne«j"  |R(e")(e'J2 
3  he  Qe 


(53) 


Note,  this  definition  is  slightly  different  than  that  used  by  other 
authors  since  the  radiation  field  term  is  not  approximated  into  Eq. 
(53)  as  was  done  previously.4^  In  our  formulation  it  is  explicitly 
retained.  Next,  we  use  the  approximate  separation  of  energy  and  the 
product  approximation  for  the  partition  function  and  transition 
moment  squared, 

Qtot  =  Qe  Qrol 

X  |R(e"J"M")(e'J'M'J2  =  |R(e")(e')P  |R(J")(J'J2 

M",M' 


Efj  =  Ee  +  Ej 
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and 


<%")(£')  =  W(£"J")(£-J') 


(54) 


Inserting  these  expressions  into  Eq.  (52)  yields 


>(e"J")(eT) 


-  8ji3  IQ-36  d£"  e~Ee''/kT 


3hc 


de"  Qe 


0)(e")(£') 


|R(e")(e')P 


dj"  dsytn  e'E],,/kT  tO(e"J")(£'J') 

dj”  dSym  Qrot  )(£) 


[rad]  |R(j”)(j  )|2 


(55) 


Writing  drotsd,jDSym  gives  the  following  expression  for  the  reduced 
line  intensity 


S(e"j")(eT) 


dr„l  e-EJ"/kT  tO(E'T)(ET) 

drot  Qrol  “(£")(£') 


[rad]  |R(j")(j-)|2 


(56) 


We  now  go  back  to  Sf\  Eq.  (53),  which  contained  a  factor 
|R(e,,)(e')l2-  This  must  be  expressed  in  terms  of  Using  Eq.(44) 

we  have 


^ 


I  A 


drivr  j j» 


eTM'— 


=  -J_  I  64^  '°'U  |R<e-J-M')teTM'f  -(57) 

dj’M'  j»tj'  3h 


We  now  assume  that  C0(e,O(e,)®a)(e,'J"M,')(e,J’M’)>  and  we  can  write 
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Ae’— >e" 


1  64n4 10'36 


dPM.  3h 


^(e-Xe)  X  |R(E"J"M“)(e7M’)|2 
J",J' 


M",M' 


=  64SM^_(03c.,vcJr 
3  h 


’(£•')(£)  lK(e"Xe')l 


(58) 


where  we  have  used  (e.g.  see  Penner45  top  of  page  136) 


IVxe'J2  =  ” —  X  lR(e"J"M")(eTM'f 
dj'M'  j-j- 


(59) 


Working  with  this  we  have, 

— —  X  |r(£"J"M"xe'J'M')|2  = — X  |R(£"xd2  |R(rxt'J2  •  (60) 

dj'M’  j"?j'  dj’M'  j'\ j’ 

MMVT 

however  ^  |R(J")(J')|2  =  dr  thus 

,r 

k-xe-J^-J-lRrxrfldr  .  (61) 

dj'M'  j- 

We  finally  arrive  at  the  result 

A(e"Xe')  =  tU(£"Xe)  “ — IR(e"X£')P  X  '  (62) 

3  h  dj'M’  j- 

Now,  we  turn  to  evaluating  the  integrated  band  intensity  according 
to 
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(63) 


Se'  =  X  S(e"J")(eT) 

J",J' 


=  X  —  10  36  0)(£"J'')(e'J')  frac,l  7“  X  |R(e"J''M")(e’J'M’)|2  -(64) 

3hc  dj-  M.,  M, 


j",r 


Rearranging  gives 

sl  =  8ftM.Q±i.  0)(e..)(e.)  X  de'T'  e^£  J  -  -J-  X  |R (e"J“M ")(£' J'M 'J2 

Q  dj- M..  M. 


(65) 


3hc 


j",j’ 


where  we  have  assumed  to  (e")(e')“  ©(e"J"M")(e,J,M,)»  [rad]~l, 
de"J"=de"dj".  Now  using  the  product  approximation  for  the  partition 
sum,  energy  additivity,  and  the  separation  of  the  transition  moment 
squared,  X  Ir(e"J"M")(etm-)|2  =  iR^xe#  |R(rxnP  ,  we  arrive  at 


?; - ^4^ ^ de"  e^E-—  m  i  -(66> 

3  he  Qe  J" ,}'  Vrot 

Next,  we  use  the  sum  rule  X  lR(J  ")(J’j2  =  dJ"  and  X  d->"  e  Ej  =  Qrot 

j,  J" 

and  we  see  the  double  sum  on  J",  J'  gives  Qrot/Qrot=l-  Thus, 

7£"  If)'36  dE"  e*Ee"/kT  Id  |2 

Se.  =  ZK  M—  CO(e")(£')  -S— - |R(e"xd 

3  he  Qe 


(67) 


Our  last  step  is  to  relate  the  integrated  band  intensity  to  the  reduced 
Einstein  A  coefficient.  Starting  from  Eq.  (57),  assuming  that 
<fl(e")(e')~C0(e"'J  "M”)(e'J'M’)’  and  X  |r(e"J"M")(eTM')|2  =  |r(e'')(e,)|2  |R(r)(J')P 


M".M' 


we  can  arrive  at 
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(68) 


A£'_>e"  -  f  CO^")(e')  l^(e")(e')|2  ““ 


3  h 


dj’M 


X  Ir(J"xj'J2 

r,J’ 


The  sum  over  J'  gives  the  sum  rule  mentioned  above, 

X  |R(J")(J')|2  =  dr  »  a°d  t*ie  sum  over  ^  "  g*ves  X  dr  =  X  (2J"+1)  which 

j'  r  j" 

is  just  the  number  of  states  in  the  lower  electronic-vibrational  state. 
The  interpretation  of  1/dj'M'  comes  from  the  meaning  of  Eq.  (44)  for 
the  reduced  quantity.  The  reduced  Einstein  A  coefficient  is  given  by 
a  sum  over  all  lower  states  of  the  Einstein  A  coefficient  divided  by 
the  total  number  of  upper  states.  Thus,  dj’M'  is  the  total  number  of 
states  in  the  electronic  vibrational  band 


dj'M'  =  Xd,'  =  X(2J'+l)  .  (69) 

r  r 


Recall,  that  for  considering  the  separation  of  the  electronic- 
vibrational  part  of  the  intensity  from  the  rotational  part,  we  assumed 
the  rotational  states  are  degenerate  in  comparison  to  the  electronic- 
vibrational  states.  Finally,  we  arrive  at 


A 


£* — >£" 


dj"M" 

dj'M’ 


(70) 


We  note  that 


X(2J"+1) 

dj"Mi  _  _£. - =  ratio  of  states  (71) 

dj-M-  £(2J'+1) 
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which  for  X^Z'g/a^Ag  is  3/2. 

The  relationship  between  the  line  intensity  and  the  reduced 
Einstein  A  coefficient  is 


S(e'T)(eT)  -  — 


dP"  dr  e'Ee"J"/kT  co(£-j-)(eT) 


87tC 


3 

“(£")(£•) 


[rad]  -J- 
dj' 


|2  dj'M' 
dj"M" 


ve'-*e 


"(72) 


and  the  relationship  between  the  integrated  band  intensity  and  the 
reduced  Einstein  A  coefficient  is 


1  dp"  e-Ee"/kT  1  _djTML  A£  _^E" 
8ftC  Qe  C0(E")(£')  dj'M" 


(73) 


7.0  HITEMP 

There  is  a  need  to  add  to  the  HITEMP  database  the  correct 
halfwidths  of  water  vapor.  This  can  be  accomplished  by  using  a 
group  of  programs  and  direct  access  files,  however,  this  must  be  set¬ 
up  for  the  particular  platform  that  the  program  will  be  run  on.  This 
is  because  of  the  nature  of  the  direct  access  files.  Some  time  ago  the 
programs  were  put  on  the  Phillips  Laboratory  VAX  where 
modifications  were  made  to  run  in  double  precision  and  other  VAX 
specific  features  were  implemented.  The  addition  of  the  halfwidths 
to  the  high  temperature  water  files  appeared  to  have  some  problems. 
Some  of  the  air-broadened  widths  were  zero,  the  extrapolated  values 
were  not  correct,  some  self-broadened  values  were  not  correct,  and 
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the  reference  and  error  codes  were  not  used.  To  remedy  this,  I  took 
the  necessary  codes  and  installed  the  ADHW-H20.F0R  program  on  a 
PC  at  Ontar  Corp.  This  meant  that  the  direct  access  files  needed  to  be 
created  on  this  computer  from  the  initial  files,  etc.  This  was  done 
and  the  program  was  tested  and  is  working  fine. 

Also  for  HITEMP,  it  was  found  that  some  of  the  transition 
assignments  are  incorrect.  Since  this  is  the  only  true  unique 

identifier  for  the  transitions,  this  must  be  corrected.  The  list  of 
corrections  was  supplied  to  us  by  Dr.  John  Selby.169  A  program  was 
written  to  check  this  list  automatically  and  search  for  duplicates. 

To  correct  the  transition  assignments  for  H2O  on  HITEMP  which 
are  incorrect,  Dr.  R.  Wattson  has  gone  through  the  list  supplied  from 
Dr.  John  Selby  with  reference  to  his  DND  programs  and  has  produced 
a  corrected  version  of  the  list.170  The  new  list  was  run  through  the 
programs  which  check  the  list  for  duplicates.  This  is  extremely 
important  because  of  the  possibility  of  dual  (or  more)  changes.  The 
danger  here  is  that  since  the  identification  is  the  label  of  the 

transition,  once  it  is  changed  it  is  changed  forever.  We  do  not  want 
to  change  an  assignment  to  something  else  and  then  find  that  what 
we  changed  it  to  will  then  change.  Our  program  found  no  duplicates. 
We  then  wrote  a  new  program  that  reads  the  assignment  change  list 
and  stores  the  information  as  a  function  of  the  quantum  assignments, 
then  reads  the  HITEMP  water  database,  finds  if  the  assignment  read 
need  to  be  changed  and  then  makes  the  change  while  updating  the 
references,  etc.  This  program  was  tested  on  a  small  subset  of  the 

data  and  by  me  and  independently  by  Dr.  Wattson  and  it  appears  to 
work  fine.  The  final  program  was  delivered  to  Ontar  where  it  was 
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run  on  the  data  and  then  destroyed  (see  comments  above).  The 
error  flags  and  reference  codes  were  also  set  by  this  program.  After 
this,  the  new  files  with  the  correct  quantum  identifications  were  run 
through  the  ADHWH20.F0R  program  to  add  the  air-  and  self- 
broadening  halfwidths,  and  the  temperature  exponent  for  air¬ 
broadening.  These  operations  completed  the  work  on  water  vapor 
for  the  HITEMP  database. 
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